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Precise Measurements With Bingham Viscometers and 
Cannon Master Viscometers 


J. F. Swindells, R. C. Hardy, and R. L. Cottington 


\ critical study has been made of the techniques used at the National Bureau of Stand- 
ards with Bingham viscometers and Cannon kinematic viscometers. All corrections appli- 
cable to measurements with these instruments were critically examined. Instruments of 
each type were calibrated using the viscosity of water at 20° C as the primary viscosity 
standard. The viscometers were used to determine the viscosities of four hydrocarbon 


liquids in the range 0.4 to 40 centipoises 


With each liquid, the values obtained in the 


two types of viscometers were in agreement by 0.05 percent or better, indicating that no 
gross error was involved in the use of either instrument. It is considered, however. that 
the inherent relative simplicity of operation of the kinematic viscometer makes it a pref- 
erable instrument for this type of measurement. 


1. Introduction 


As the result of a recent determination [1] ', the 
National Bureau of Standards on July 1, 1953, 
adopted the value of 1.002 centipoises (cp) for the 
absolute viscosity of water at 20° C as the primary 
standard for viscosity determinations. The Ameri- 
can Society for Testing Materials, The National 
Physical Laboratory in England, and the Physika- 
lisch Technischen Bundesanstalt in Germany have 
concurred in this action. Previous to this, the 
values of the secondary standards of viscosity issued 
by the Bureau were based upon 1.005 cp for the 
viscosity of water at 20° C. In connection with the 
reevaluation of the secondary standards on the basis 
of the new value for water, a comparative study has 
been made of the use of two types of viscometers for 
relating the viscosities of other liquids to that of 
water as a primary standard. Bingham viscometers 
and Cannon viscometers were used, and comparisons 
were made of the viscosities of four hydrocarbon 
liquids in the range 0.4 to 40 cp and of the viscosity 
of each liquid as determined in each type of instru- 
ment. This paper presents in some detail the 
techniques used in making these determinations to 
describe the methods employed in evaluating the Na- 
tional Bureau of Standards standard viscosity sam- 
ples and to call attention to the magnitude of certain 
corrections often neglected in viscometry. The ex- 
tension of these techniques to the calibration of 
viscometers with larger capillaries suitable for the 
measurement of the viscosities of more viscous 
liquids is relatively simple and involves the same 
methods as are covered here. 


2. The Bingham Viscometer 


2.1. General 


A short treatment of the use of the Bingham 
viscometer, shown diagrammatically in figure 1, has 
been given previously [2], and it will be assumed that 


| Figures in brackets indicate the literature references at the end of this paper 


the reader has some familiarity with the instrument. 
With the usual procedure, pressure is applied to the 
right limb of the viscometer, and the rate of flow is 
determined by measuring the time required for the 
meniscus in the left limb to pass from the fiducial 
mark d to mark c. To avoid the necessity for drain- 
age corrections, each determination is made with the 
bulb A initially dry. This is accomplished by intro- 
ducing the sample into bulb B with a special pipet, 
sufficient liquid being added to fill the viscometer 
between the marks d and g at the test temperature. 
With the exception of certain calibration runs with 
water, the viscometer is cleaned and a fresh sample 
is introduced for each measurement. By this pro 
cedure the volume of flow is kept constant for each 
instrument regardless of the viscosity or rate of flow 

Pressure is applied to the liquid in the viscometer 
by air supplied from a tank having a capacity (about 
60,000 em*) sufficiently large that the increase in vol- 
ume (about 4 cm*) in the pressure system during the 
flow of liquid from bulb B causes no significant re- 
duction in pressure. The tank is thermally insu- 
lated to prevent rapid changes in its temperature 
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Figure 1 Bingham viscometer, 
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with fluctuations in the ambient temperature. In 
general, pressures above 150 mm Hg are read with a 
mercury manometer, a water manometer being used 
for lower pressures; differences in liquid levels in the 
manometer up to 600 mm are read with a cathetom- 
eter and greater differences are read by using a cal- 
ibrated steel tape and reading The 
heights were read to about 0.02 mm with the cathe- 
tometer and 0.1 mm with the steel tape 


telescopes. 


Times of flow are measured using a stopelock oper- 
ated from a constant-frequency source of power 

All of the flow measurements reported here were 
made with the viscometers immersed in a well-scirred, 
30-gal oil bath whose temperature was controlled by 
supplying a heater from the output of a manually 
adjusted variable transformer. When operating be- 
low ambient temperature, heat was taken from the 
bath through cooling coils in which cold water was 
circulated at a constant rate. ‘Temperatures were 
measured with a resistance thermometer, and their 
recorded values are believed to be accurate to 
-0.002° C, 


2.2. Calibration 


The usual form of the modified Poiseuille equation 
for the calculation of viscosity by the capillary 
method is 


rr'Pt mpV 


—_ ; _ , (1) 
BV i+ nr) Sxil+nnt 


where 


radius of capillary 

mean effective pressure drop through the 
capillary 

volume between fiducial marks 

time for volume V to flow 

length of capillary 

coefficients associated with the flow at the 
ends of the capillary 

the absolute viscosity and density of the 
liquid whose viscosity is to be deter- 
mined; 


As over a considerable range of rates of flow, m and n 
can be taken as constant for capillaries having square- 
cut ends [3], certain of the quantities in eq (1) are 
usually grouped to give two constants for the instru- 
ment. Equation (1) then takes the form 


2=CPt—C’ pit, 
where 


. mV 
and C’ =. . 
Sx(l+nr) 


> =ra 
~ 8V(l+nr) 


With viscometers having capillary bores small enough 
for calibration with water, these two constants are 
commonly evaluated by measuring the times of flow 
of water at 20° C for various applied pressures. 
The product Pt is plotted against 1/t, and C and C’ 


are evaluated from the resulting straight line by in 
troducing known values for 7 and p. For the tw 
water-calibrated viscometers used in the work 1 
ported here, however, it was found that such plot 
did not yield straight lines over the complete rang 
of rates of flow in which the instruments were to bh 
used, indicating that m and n were both varying, an 
consequently C and C”’ could not be treated as cor 
stants. For this reason and for simplicity, t) 
purely empirical formula, 


g=fFt, 


was used with the viscometers. In this equation, 

is a multiplying factor that varies with the condi 
tions of flow. It is convenient to assume that th 
conditions peculiar to a given rate of flow are char 
acterized by the corresponding value of the Reynolds 
number as calculated for the flow in the capillary 
(R=2Vp/xrnt), and therefore f will have a definite 
value for each value of the Reynolds number. The 
calibration of these instruments, then, consists of 
determining the value of f as a function of the Reyn 
olds number. 

Although the Bingham viscometer is designed to 
minimize the effect of hydrostatic head in the 
instrument, actually small head corrections to an 
applied external pressure, po, are necessary to give an 
exact value of P. The first correction, which will be 
called the level head, arises from the fact that, despite 
careful construction, the two bulbs (A and B in fig. 1 
will not be geometrically identical and at the same 
level. This condition results in a residual hydro- 
static head that is of the same magnitude but 
opposite sign for flow right to left (R-L) and left to 
right (L-R) in the viscometer. The correction is 
largely independent of the magnitude of the applied 
external pressure but is proportional to the density 
of the liquid. The level head was evaluated in this 
work by making flow measurements in the L-R 
direction and then repeating in the R-L direction, 
using approximately the same applied pressure. The 
runs were made at the lowest applied pressure that 
could be accurately measured, and under conditions 
chosen to insure that the volume of flow was the 
same in each direction. Under these conditions, the 
difference in the time of flow for the two directions 
was attributed to a corresponding difference in 
pressure. Then, knowing the applied pressures, the 
magnitude of the correction was calculated. 

A second head correction, which becomes very 
significant at low values of the applied pressure, is 
the commonly termed logarithmic head correction 
If we neglect for the moment the level head discussed 
above, the initial pressure drop through the capillary 
is Po+ax, where z is the initial hydrostatic head of 
liquid, and the final pressure is po>—z, but the mean 
effective pressure is not exactly po, the arithmetic 
mean of the two, since the higher pressures during 
the first part of the run cause more rapid flow and 
are therefore effective for a proportionately shorter 





For bulbs of regular shape, however, the dry-bulb runs made with the same applied pressure 
etive pressure is readily calculable. For this On this basis, all the wet-bulb runs were corrected to 
k, the caleulations were made by using Barr’s  dry-bulb conditions. 
ation [4]. This equation is written P(1+ "7/10 The calibration data obtained with water at 20° C 

,5 : )=Po, fapproximately, where y for viscometers | and 20 are given in appendices 7.1 

Barr states that this equation is derived “for and 7.2. All runs were made in the R-L direction. 
ulb that has the form of a pair of opposed cones, except for the few special runs made for the deter- 
of a height equal to the radius of the common mination of the level head corrections. From these 

discharging into air, or for such a bulb dis- data the plots of Pt versus R given in figures 2 and 3 
rging symmetrically into a similar one, as in were constructed. Each plotted point represents 
ngham’s viscometer’. The correctness of this the mean of at least two runs made at approximately 
yression was verified in an independent derivation, the same pressure. For each viscometer, it is seen 
because the radius of the common base dropped that the plot is not linear at certain values of R, 
in the derivation, it is concluded that the ex- as was mentioned previously. Since the data are 
ion holds for cones of any height-to-base ratio. not adequate to position accurately the curves at 
of importance to note in connection with runs — the lowest Reynolds numbers, the lines were drawn 
low pressures, that the logarithmic head is a horizontal in this region on the basis of theoretical 
ction of the ratio of the initial and final pressures, considerations [5]. From these curves values of 

d therefore it increases in absolute magnitude as Pt corresponding to selected values of R were ob- 

applied pressure is reduced. The sign of the — tained and, taking 7=0.01002 poise, were substituted 
rection is always negative, no matter which 
ection of flow is being used in the viscometer, and 
erefore it does not affect the determination of the 
vel head, as described above 

The design of the Bingham viscometer is such as 
to minimize the effect of surface tension upon P 
Calculations of the net contribution of surface tension 
to the mean effective pressure indicate a maximum 
correction of about 0.001 mm Hg, which was ne- 
vleeted in this work. 

\ more complete discussion of the above correc- 
tions is given by Barr [3, chap. 3] 

The three Bingham viscometers used are identified 
in table 1. Of these, viscometers 1 and 20 were 
calibrated with water, whereas number 5 was cali- 
brated subsequently by using two hydrocarbon oils 
whose viscosities were determined in viscometers 1 
and 20. In making the calibration runs with water 
an exception was made to the usual procedure of 
flowing into a dry bulb, most of the runs being made = Ficure 2. Pt as a function of Reynolds number for Bingham 
with the fiducial bulb wet. Following each R-L viscometer number 1. 
run, the liquid war drawn back at such a rate as to 
require about 80 see with viscometer No. 1, and 120 
sec with No. 20 for the fiducial bulb to empty. After 
allowing about 18 min for further drainage from the 
walls of the bulb, another R-L run was started. In 
this manner a reproducible volume of flow was 
obtained without the necessity of drying the instru- 
ment between runs. The difference between the 
volume of flow under this wet-bulb condition and the 
true volume of the fiducial bulb was obtained from 
the differences in the product Pt between wet and 


ee et 








SIRI OTT, FOS SFT tt ee 
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Ficure 3 Pt as a function of Reynolds number for Bingham 
viscometer number 20 
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in eq (3) to calculate corresponding values of /f. 
The values of f were then plotted as a function of R, 
vielding the calibration curves for viscometers | and 
20 shown in figures 4 and 5, respectively. 

To calculate the unknown viscosity of a liquid 
from flow data obtained in one of these viscometers, 
it is necessary to choose a value f from the calibration 
curve corresponding to a tentative value of the 
Reynolds number based on an estimate of the vis- 
cosity of the unknown. If the viscosity value 
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calculated by substituting this value of f, togeth 
with the measured values of P and ¢, in eq (3) diffe: 
materially from the first estimate of the viscosit 
a second calculation is required. 

Viscometer No. 5 was calibrated with oils A an 
B (to be identified later), using values of their vi 
cosities as determined in viscometers 1 and 2 
These calibration data are given in appendix 7 
As all calibrating runs and viscosity measuremen 
made in viscometer 5 were at Reynolds numbe 
less than 10, it was assumed that f in eq (3) wa 
constant for this viscometer under these conditions 
The level head was measured and the logarithmi: 
head calculated for this viscometer as outlined 
previously. 


2.3. Results with Bingham Viscometers 
The viscosities of normal heptane and_ thre 


additive-free hydrocarbon oils were measured at 
20° C. The mean results are given in table 2, and 


Prarie 2 Results at 20° C with Bingham viscometers 
Viseo » 
oi enetes / ” 7. mean 
p j 

#- He ptare l “YY 0, 4112 | 

Do m) 1M.3 gio 6of)| OM AEE 
\ 1 W.2and 24 l a) ) - 
4 » 18.1 1.9251 jf |. 9250 
B 2.4 7.09 | - 10 
8 » 1.5 7.610 j ‘61 
( 0.2 and 0 42.82 42.82 


data from individual runs in appendix 7.4. Normal 
heptane was included to furnish comparative results 
with a liquid whose viscosity is less than that of 
water. The heptane was not of highest purity, but 
was taken from a lot which met specifications for a 
primary reference fuel for the determination of octane 
number [6] and was at least 99.5 percent pure 
The viscosities of n-heptane and oils A and B were 
determined in viscometers 1 and 20, employing the 
normal procedure of running in the R-L direction 
with the fiducial bulb initially dry. Oil A was run 
at as high and as low Reynolds numbers as was 
practical. Higher applied pressures were avoided 
because experience has shown that oils as light as A 
dissolve air rather rapidly under these conditions, 
with a consequent lowering of their viscosities 
Conditions for the other runs were chosen to strike 
a balance between optimum conditions for accurate 
measurement of pressure and of time. In calculat- 
ing the viscosities, small corrections were made to 
account for the change of viscosity of the oils with 
pressure. Observed viscosities were reduced to a 
pressure of 1 atm, as described elsewhere [2], by 
adding a second term to eq (3), which then becomes 


n={Pt—FPn, (4 














e F isa factor representing the fractional change 
, for unit change in pressure. Good agreement 
=; found between results obtained in the two 
ometers 1 and 20, the largest discrepancy being 
5 percent with n-heptane. In addition, results 
| ained with oil A in viscometer 1 at R=28.2 and 
24 were in equally agreement ap- 
ndix 7.4), which tends to confirm the shape of the 
bration curve at the lower range of Reynolds 
mbers (fig. 4). 
Runs made with oil C in viscometer 5 at high and 
low pressures were In good agreement (see appen- 
7.4), indicating that the corrections applied for 
ange of viscosity with pressure are not in serious 
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or 
3. The Cannon Viscometer 
3.1. General 


In outline, the techniques used at the National 
j Mureau of Standards with the Cannon viscometers, 
hown in figure 6, are substantially as described by 

innon {7], the principal differences consisting of a 

ore rigorous treatment of the various corrections 

to be applied to the observed data. With instru- 

ents of the kinematic type the value of P applicable 
neglecting surface tension effects, is 








to eq { l 
ven by 
P hg! p Pa), (5) 


which A is the mean effective head, p is the liquid 
density, and p, is the mean density of the air column 
i the left arm of the viscometer, as shown in figure 6. 
For application to the Master viscometers, eq (1) 
then becomes 


wry ) mV 
: ait h(1 Faye - (6) 
p SVil+nr) p Sa(l+-nr)t 
If we let 
4 
nrregd Pa ° 
> s h ] - K 
8 V(/+-nr) ( p ) 
nr 
m V > 
B, 
: Sar(l+-nr) 
eq (6) becomes 
— = 
Kt——- (7) 
p t 


itis apparent that A and B will not be constant for 
all test conditions, but their values will reflect any 
variations in m and n with Reynolds number. As 
s done with the Bingham viscometers, these changes 
m and n are implicitly determined experimentally 

is a part of the viscometer calibration. _ In addition, 
the value of K changes with any change in A or 
p/p), and a small correction for variation in the 
ilue of g may be required if the instrument is used 


JP PER CAPILLARY 
SS 




















ONS IN Cm 





f Omens 


Vaster 


viscomete 


Cannon 


Ficure 6 


at a location other than where it is calibrated. If 
we let AK, be the determined value of K under the 
particular hydrostatic head at calibration, we can 
find from it the value of AK under other test condi- 
tions by calculating changes in Ah and (1—p,/p). 
Since a change in (1—p,/p) is equivalent to a change 
in A, it is grouped with the other factors which 
affect A. For any test condition, eq (7) may then 
be written 


3 = K. ( i+% i=, (8) 


p 


in which Ah is the difference in effective head between 
test and calibrating conditions, and is made up of 
the sum of Ah,, Ahs, ete., arising from the several 
factors involved. In practice it is convenient to 
treat K as a constant equal to K,, and to apply the 
corrections to the time of flow t. The viscosity is 
calculated by means of the equation 


B 


Ki.— "> (9) 


in which t.=t(1 +-Ah/A). 

Of the factors affecting A, the filling volume will 
be considered first. To introduce a reproducible 
volume of liquid into the viscometer, the viscometer 
and the liquid must both be at some standard tem- 
perature (20° C for this work) each time the vis- 
cometer is charged. To obviate the inconvenience 
of this procedure, fillings are made at ambient tem- 
perature, which is observed, and an appropriate 
correction, Ah,, is then calculated to account for 
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any difference between the filling at the observed 
temperature and a standard filling at a temperature 
of 20° C. The calculation of this correction is 
similar to that of Aho, which follows 

A second factor affecting A, which gives rise to a 
correction Ah», is the change in effective head with 
differences between the calibration temperature and 
that of subsequent runs. The calculation of this 
correction is based on the cubical expansion of the 
viscometer and test liquid, together with the mean 
diameter of the working part of reservoir C, figure 
6, |S] For the most precise work an accurate 
know ledge of the internal diameter of C is necessary. 
This measurement was made by plugging the lower 
end of the capillary with wax and then filling the 
bottom part of the viscometer with mercury to a 
low level in C. Weighed increments of mercury 
were added, and the resulting increases in level 
in C were measured with a depth gage, from which 
data a mean internal diameter was calculated 

From differences in the factor (/—p,/p) between 
calibrating and operating conditions a Ah, is calcu- 
lated. The density of the air in the viscometer is a 
function of the temperature and amount of water 
vapor present, neither of which is known accurately, 
but no serious errors are introduced by assuming the 
air to be 50-percent saturated and at the test tem- 
perature. 

The factor most difficult to evaluate is the varia- 
tion of the mean head with surface tension. In the 
work reported here, viscosity measurements of oils 
made in water-calibrated viscometers required the 
most careful estimates of the effects of surface ten- 
sion. In the case of water, the contribution of 
surface effects to the effective head was determined 
experimentally and also calculated. The experi- 
mental determination was made with a glass bulb 
blown as nearly as possible to the size and shape of 
the fiducial bulb in the viscometer. The bulb was 
connected through a U-shaped tube to a eylindrical 
glass tube about | cm in diameter and held vertically 
close to the side of the bulb. The bulb and tube 
were then mounted in a small bath of oil having close 
to the same index of refraction as the bulb. Water 
was introduced into the tube until the level stood in 
the neck at the bottom of the bulb. The difference 
between the liquid levels in bulb and tube was then 
measured with a micrometer microscope mounted 
for vertical measurements. From this difference 
and a calculated value for the capillary rise in the 
cylindrical tube, the rise at this particular level in 
the bulb was obtained. The calculation of the rise 
in the cylindrical tube was made by using Sugden’s 
table {9}. A known volume of water was then added 
and the measurements repeated. The portion of the 
added volume going into the bulb was the difference 
between the added volume and the volume increase 
in the tube as calculated from the change in level in 
the tube. This process was repeated, taking about 
20 increments of volume to fill the bulb. From 
these data and a measured value of the head in the 
viscometer at a known level in the bulb, a value of 


the mean head in an interval, A,, was estimated { 


each of the 20 intervals. Then, knowing the volun 
AV, and the capillary rise, A,, for each inter, 
terms of the form AVA,/h; were calculated for ea 
interval. The quotient 
wi AV 
— oi 


Al 


sr hy 


h, 


then gave a time-weighted mean value for the cap 
larv rise in the bulb as a whole. The rise in thy 
evlindrical lower reservoir was calculated by usi: 
Sugden’s table and subtracted from the mean valu 
for the bulb to get a net effect for water in th 
viscometer. The same techniques were 
determine a net value of the surface tension for each 
of the oils tested. From the difference between 
water and oil, a correction, Ah, was calculated for 
application to data obtained with oil. 

The net surface tension effect for water in the 
viscometer was also obtained by calculation alone. 
For this purpose, a vertical section through the 
fiducial bulb was plotted on a large scale as accurately 
as possible. The bulb was then divided into 13 
sections, and a value of AVA,/h, was calculated for 
each. With the aid of Bashforth and Adams’ tables 
110], each value of AV’ was calculated by using the 
geometry of the bulb as plotted and the calculated 
change in meniscus volume between the top and 
bottom of the section. The tables were also used 
in calculating h. for each section. This method has 
been given in more detail by Barr [11]. The mean 
value for the bulb was obtained weighting each 
section with respect to time, as was done in the 
experimental method. The calculations yielded a 
value within 0.01 percent of the experimental result, 
which lent confidence to both methods. Because 
of limited range of Bashforth and Adams’ tables, 
the calculation could not be made for the oils. 

Thus far in the treatment of the factors affecting 
K in eq (7), it has been assumed that the volume of 
flow, V, is constant under all conditions. For 
application to the work reported here, the validity 
of this assumption was carefully tested. Bulbs 
similar in size and shape to the fiducial bulbs 
in the Cannon viscometers were blown with capil- 
lary stems about 1 cm long above and _ below 
the bulb. To detect differences in the volume of 
flow of water and the several oils from a bulb, 
a capillary was attached below the bulb with neo- 
prene tubing so that the pulb could be quickly dis- 
connected. The capillary bore was selected to give 
about the same time of flow for a particular liquid 
as was found in the Cannon viscometer. By weigh- 
ing the test bulb empty and again after a liquid had 
been run from it at the rate controlled by the 
resistance of the capillary, the volume of liquid 
remaining on the walls of the bulb was obtained 
From the difference between results obtained with 
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could be calculated. Only one instance of a 
rection as large as 0.01 percent was found, and 

s so-called drainage correction was applied. 

The viscometer is filled by holding it in an inverted 

osition with the end of the upper capillary immersed 

the liquid under test. Liquid is drawn into the 
cometer by applying suction to the other arm 
il the level reaches a position slightly beyond the 
ark F (fig. 6). Applying just enough suction to 

ld the liquid level at approximately this position, 

upper capillary is lifted from the liquid and the 
wrator’s finger is placed over its end. With the 
strument still inverted, excess liquid is carefully 
ed out of the capillary past the finger tip until 
liquid level reaches F. The viscometer is then 
turned to its normal upright position, and the finger 
removed with a gentle wiping action. In this man- 
er a reproducible filling volume is introduced. 
\fter filling the viscometer it is mounted in its bath, 
nd sufficient time is allowed for the fiducial bulb 
is to empty and for the sample to come to the bath 

emperature before starting a test. To make a 
lctermination, the liquid is pushed up under 17- 
mm-Hg air pressure until the upper meniscus stands 
ibout 5 mm above the mark at the top of bulb B 
lhe air pressure is then released, and the time for 
bulb B to empty is observed. With relatively low 
viscosity oils, raising the liquid with vacuum at the 
start of a run sometimes results in a time of flow 
slightly higher than when air pressure is used. Such 
differences are attributed to the removal of aur from 
ihe sample under reduced pressure. It 1s therefore 
onsidered better technique to use air pressure as 
any change in the amount of air dissolved in the oil 
will start at the surface farthest removed from the 
liquid passing through the capillary. 

When the liquid is pushed up at the start of a run, 
the meniscus level in reservoir C is lowered, leaving 
some liquid behind on the walls. The amount of 
liquid left on the walls influences the level in C and 
hence the value of A. By always using the same 
pressure in pushing the liquid up, the level in C 
is lowered at a rate inversely proportional to the 
liquid viscosity, which leaves about the same amount 
of liquid on the walls for each test. Calculations 
indicate that this procedure eliminates significant 
variations in h. 


3.2. Calibration and Results 


The calibration of the Cannon viscometers con- 
sists in evaluating K, and B in eq (9). To illustrate 
the method, eq (9) is rewritten 


ae 
- =K.— 3) (10) 
pt. | 


n which @ is substituted for ¢4 without detectable 
Liquids having different kinematic viscosities 
n/p) are run in the viscometer, and values of 7/(pf,) 
ire plotted against 1/f. Values of A, and B are 


alibrating liquid and a test liquid, a correction to 
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then obtained from the resulting curve. As un- 
published results in this and other laboratories have 
shown that B may not remain constant over the 
useful range of Reynolds numbers for this type of 
viscometer, the calibration must embrace the range 
in which the instrument is to be used. In practice, 
therefore, the calibration curve consists of a plot of 
n/(pt.) asa function of 1/@ or 2, and for a particular 
viscosity determination values of A, and B are 
obtained from the curve at the appropriate value of 
l f? or R. 

The four Cannon viscometers used in this work are 
identified in table 3. Of these, M25-1 and M25-2 
were calibrated with water, whereas \1104 and M105 
were calibrated subsequently with two hydrocarbon 
oils whose viscosities were determined in M25-—-1 and 
M25-2. The calibration data obtained with M25-1 
and M25-2 are given in appendix 7.5, and from 
these data the calibration curves in figure 7 were 
plotted. 

The curves were arbitrarily drawn through the 
values of »/(pf,) at Reynolds numbers of about 15, 
which were obtained with water at 20° C, because a 
value for the viscosity of water at 20° C was taken 
as the primary standard for calibration. The points 
at higher Reynolds numbers were obtained with 
water at 40°, 60°, and 80° C. The values used for 
the properties of water at these temperatures are 
given in appendix 7.6. As A, refers to the conditions 
existent with water at 20° C, times of flow obtained 
at the higher temperatures were corrected by appro- 
priate values of the factor (1+ Ah/h) in eq (8). 

It is evident that the value of B is zero (m—0) for 
these two viscometers up to a Reynolds number of 
at least 110. Because they were used at values of 
R never exceeding 50, the determinations made in 
these instruments were calculated by means of the 
equation 
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TaBLe 3 Essential dimensions of Cannon viscometers 
> Filling Radius of 
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M252 4.5 Obs $13 “6.3 t 150 
Mim “4.3 24 i wif 6s 1.53 
Mit “4 rou MM wf ; 1 


Viseometers M104 and M105 were calibrated with 
oils A and B, using values of their viscosities as deter- 
mined in M25-1 and M25-2. The calibration data 
are given in appendix 7.7 These viscometers were 
calibrated and used at R<10, and hence it was 
assumed that B=0 and eq (11) was applicable. 
For consistency, the values of A, for these viseom- 
eters were also calculated for the conditions that 
would exist with water, even though water was never 
run in these viscometers. The values of the viscos- 
ites of n-heptane and the three mineral oils, as deter- 
mined in the Cannon viscometers, are given in table 
4. Detailed calculations are given in appendix 7.8 
In table 4, kinematic viscosities are converted to 
absolute units for comparison with the measurements 
made in the Bingham viscometers. The results 
obtained independently in the two types of viscom- 
eters are in good agreement, the largest discrepancy 
being 1 part in 2,000 in the case of n-heptane. As 
the techniques employed and the nature of the cor- 
rections applied are quite different with the two types 
of instruments, it is not likely that similar errors are 
involved in both methods. Hence, there is indica- 
tion that the treatments of the data obtained in both 
types of viscometers are free of gross error. 


TABLE 4 Results at 20° C with Cannon viscometers compared 
with Bingham results 
Viscom n. by 
a) »e > ) " 
ul eter ’ we meen ‘; Bingham 
cs cs giem cp ep 
n-Heptane Ds 7 . 
= anne pone 1 0 am } 0. e00TO 0. 6841 0. 4100 0.4111 
M25 2% 
7 + i) 2 ro 2. 4043 TS110— «1. 9240 1. 9250 
4 M25 0. SK - > 
oae-2 | © eis j % Inte Cases os 7.610 
Cc M104 19. 583 . » 
‘ M105 19. 585 } 19 mt seit? 42.83 42. 82 


4. Discussion of Results 


4.1. Sources of Error in Using Bingham Viscometers 


From a thorough examination of the various fac- 
tors involved in the use of the Bingham viscometers, 
it appears that the greatest source of uncertainty 
arises from the unavoidable use of relatively low 
pressures in establishing the portion of the calibra- 
tion curve corresponding to low Reynolds numbers. 
The lower pressures are more difficult to maintain 


relatively steady and measure with a given relatj 
precision, and larger errors are introduced for giv, 
uncertainties in the values of the head correction 
This point is illustrated from the data for viscomet; 
number 1 by a comparison of run 1-24, made at 


Reynolds number of 8.5, with run 1-38 made 
R=41 (see appendix 7.1). For run 1-38, with a 
applied pressure of 73 mm Hg (100 em H,O), ¢! 
level head correction is 0.07 percent, and the logarit 
mic head correction is 0.01 percent, whereas for ru 
1-24 with an applied pressure of 15 mm Hg (21 «1 
HO), the level head correction is 0.31 percent, an 
the logarithmic head correction is 0.24 percent 
Thus in the necessary extension of the calibratioy 
curve from R=41 down to R=8.5, the sum of th 
head corrections is increased from 0.08 percent to 
0.55 percent of the applied pressure. For each of 
the three Bingham viscometers used in this wor! 
the value of the level head was determined expe: 
mentally, as described previously, the accuracy of 
the determinations being limited chiefly by the 
precision with which the pressures could be measured 
The experimental evaluation of the logarithmic head 
corrections with sufficient accuracy, however, was 
not found feasible, and hence, the corrections wer 
calculated, use being made of an equation given by 
Barr [4], by which the mean effective pressure for a 
run is calculated from the known initial and final 
pressures. This equation was developed on th 
assumption that each bulb of the viscometer had 
the shape of two cones placed base-to-base.  Al- 
though this simple geometry is not exactly realized 
in any of these viscometers, the correction is not very 
sensitive to instrument dimensions. Further confi- 
dence is lent to the calculations by the fact that good 
agreement is always found between viscosity meas- 
urements made first under conditions where the 
correction is relatively large and again when the 
correction is small. This is shown to some extent 
by the viscosity determinations made on oil A in 
viscometer 1 (see appendix 7.4). This oil was run 
first at a pressure of 20 mm Hg, with a logarithmic 
head correction of 0.08 percent, and then at a pres- 
sure of 237 mm Hg, where the correction is negligible 
The agreement between these runs is very good 
which indicates that the shape of the calibration 
curve, and the values of the head corrections applied 
are essentially correct. 

It has been shown [12] that the viscosity of water 
is not significantly affected either by change in __ the 
amount of air in solution or by change in pressure 
under the conditions of these tests. Although the 
viscosities of the hydrocarbon liquids are mor 
sensitive to these variables, the experimental con- 
ditions were chosen such that the variations in thy 
dissolved air under the different test conditions 
would not be expected to have measurable effects 
Furthermore, unpublished experiments at the Bureau 
have shown that the corrections as applied for th 
change of viscosity of these hydrocarbon liquids 
with pressure [2] are probably not in error by mor 
than 30 percent. In the tests recorded in Appendix 
7.4, the largest correction occurs with runs 5-17 am 
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s on oil C, with an applied pressure of 595 mm 

As the pressure correction for these runs is 
S percent, an error of 30 percent in the correction 
amount to only 0.02 pereent im the Viscosity 
rmination, 


Sources of Error in Using Cannon Viscometers 


he first considerations are concerned with the 

iblishment of the calibration curves given in figure 

wr the two viscometers calibrated with water. In 
work reported here these viscometers were used 

Reynolds numbers in the range 0.03 to 42, and 
nce only factors affecting this portion of each curve 
| influence the results. As the determinations in 

viscometers were to give results based on a value 

the viscosity of water at 20° C, each calibration 
rve was arbitrarily put through the point at a 
vnolds number of about 15, which represents the 
ta obtained at this temperature. At this point, 
certainties in the value of AK,, which is equal to 
of, will depend only on errors in ¢,, as the value of 
s taken without error and p is known with high 
ccuracy. With reference to appendix 7.5, it is seen 
hat for the data at 20° C, ¢, is obtained from the 
‘bserved time of flow ¢, by applying only a relatively 
small correction for filling temperature. This cor- 
rection is readily calculable and is not considered as 
a source of error. Of the factors affecting ¢, the 
effects of errors in temperature measurement, or of 
lemperature variations, are probably small. Both 
the platinum resistance thermometer and the Mueller 
(;-2 bridge used for measuring temperature were 
calibrated shortly before or after tests were made, 
so that no error larger than 0.001 deg C is expected 
i the measured temperatures. Observed tempera- 
ture variations during the tests were 0.002 deg C or 
less, with the mean temperature closer than is indi- 
cated by this spread. It is probable that the actual 
bath temperatures were within 0.002 deg of 20° C 
which is equivalent to an uncertainty of only 0.005 
percent in the viscosity of water. Calculations’ of 
the precision of the mean values of ¢ as recorded in 
appendix 7.5 indicate an uncertainty in the mean 
time of flow for viscometer 125-1 of +0.03 percent 
and for viscometer M25-—2 an uncertainty of +0.007 
percent. These limits reflect errors due to variations 
in filling the viscometers, bath temperature, timing, 
and unrecognized factors. From the above consid- 
erations, the values of K, at 2=15 for the two vis- 
cometers are probably within +0.03 percent of their 
true values. 

Having established one point on each calibration 
curve, the other important factor is the shape of the 
curves in the range R<50. From theoretical con- 
siderations [3], with increasing values of FR these 


curves should be horizontal lines until some value of 


’ is exceeded and a gradually increasingly negative 


slope begins to develop. No known theory can 
count for a positive slope at these Reynolds num- 


The calculations were made by the American Society for Testing Materials’ 
thods [13], assuming that the true value of the mean will lie within the limits 
en 99 times in 100. 
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bers. The values of »/pf. for water at 40°, 60°, and 
80° C were therefore obtained and plotted to indicate 
whether the curves could be assumed to have zero 
slope up to R=50. These values of »/pt. were com- 
puted by using the values of n/p given in appendix 
7.5 and are believed to be accurate within 0.1 pereent 
It is seen that zero slope is indicated in the range of 
Reynolds numbers covered, and unreasonably large 
errors in the values of n/pft. for water at 40°, 60°, 
and 80° C would be required to change the positions 
of the calibration curves at R<50 by as much as 
0.01 percent 

The calculations of the viscosities of n-heptane and 
oils A and B from the data obtained with these two 
water-calibrated viscometers are recorded in appen- 
dix 7.8. These data show that the times of flow for 
repeat fillings of a given viscometer with a given oil 
are in even better agreement than was the case with 
repeat fillings with water. It is estimated that the 
mean times of flow for the two fillings with each oil 
should be within +£0.01 percent of their true values, 
based upon 99-percent confidence limits. 

In appendix 7.8 the values of ¢, are obtained by 
applying at least three corrections to the observed 
times of flow. Of these, Ah,/h is calculated without 
significant error. With reference to section 3.1 it is 
seen that Ahh is dependent upon differences in 
p air/p liquid between calibrating and test conditions. 
The values of Ah,/h were calculated by assuming the 
air in the viscometers to be 50-percent saturated 
under all conditions, which was not always the case. 
No significant error is introduced by this assumption, 
however, as the difference between the values of 
Ah,/h, assuming dry air and again assuming satura- 
tion, affects ¢, by only 2 parts in 100,000. The third 
correction, Mig/h, which ts applied to account for the 
difference in head caused by the difference in surface 
tension between water and the oils, is perhaps the 
greatest source of uncertainty. As described previ- 
ously, the calculated value of the effeet of surface 
tension for water at 20° C was in good agreement 
with the experimental value, but both calculations 
and experiments were sufficiently complex as to 
make it difficult to estimate the accuracy of either. 
It is believed no error greater than +0.02 percent is 
introduced by this correction, but this can not be 
said with certainty. 

The drainage factors recorded in appendix 7.8 
based upon experiments described in section 3.1, are 
believed to be known to about 0.002 percent and are 
therefore not significantly in error. No satisfactory 
explanation has been found for the different behavior 
of oil A, as compared with n-heptane and the other 
two oils, in bulbs of similar shape and volume to the 
fiducial bulbs of viscometers M25-1 and M25-—2 (see 
fig. 6). Oil A did not show this different behavior 
in a somewhat more elongated bulb similar to the 
fiducial bulbs of viscometers M104 and M105. 

In the calibration of viscometers M104 and M105 
with oils A and B and the subsequent determinations 
of the viscosity of oil C, the sources and magnitudes 
of error are essentially the same as described above 
for the oil runs in viscometers M25-1 and M25-2. 
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4.3. Heating Effects in Capillaries 


So far in the treatment of the data it has been as- 
sumed that, once the liquid under test has reached 
the temperature of the bath in which the viscometer 
is immersed, it remains at this temperature during 
the measurement. Because heat is generated in the 
liquid as it is sheared in the capillary, it is obvious 
that this assumption cannot be strictly true. As the 
temperature will vary from point-to-point within the 
capillary, the viscosity will vary accordingly, and 
what is measured is a mean value of the viscosity cor- 
responding to some effective mean value of the tem- 
perature within the capillary. In attempting to 
calculate this mean effective temperature, one is un- 
able to postulate what conditions of temperature dis- 
tribution and heat transfer exist for a given flow con- 
dition. Hersey [14] and Hersey and Zimmer [15] 
have derived equations for calculating the mean 
effective temperature rise, first assuming conditions 
of adiabatic flow and then for conditions of thermal 
equilibrium in the liquid flowing in the capillary. 
Corrections based on these equations have been cal- 
culated for the data in both the Bingham and Can- 
non viscometers., Referring to appendices 7.4 and 
7.8 it is seen that oil C was run in a Bingham vis- 
cometer at pressures of 595 mm Hg (runs 5-17 and 
5-18) and 303 mm Hg (runs 5-15 and 5-16), and 
also in Cannon viscometers at 29 mm He. If we 
should assume the conditions defined by Hersey and 
Zimmer as incomplete adiabatic flow and calculate 
corrections, using their eq (5) 115], we get a Viscosity 
correction of +-0.70 percent for the tests in the Bing- 
ham viscometer at 595 mm Hg, and a correction of 

+0.02 percent for the tests in the Cannon viscom- 

eters. On this basis the tests at 595 mm Hg uncor- 
rected, should be about 0.7 percent lower than the 
runs at 29mm Hg. As the results recorded in table 
4 and appendix 7.4 are in good agreement at all three 
pressures without these corrections, it must be con- 
cluded that conditions of incomplete adiabatic flow 
are not approximated in these tests. Calculations 
based on Hersey’s eq (13) [13], which assumes flow 
under the conditions of thermal equilibrium, vield 
negligibly small corrections, in agreement with the 
evidence of the tests referred to above. 

It should be noted that the effects of heating in 
the capillary are not only probably negligibly small 
under the conditions of the tests reported here, but 
in addition, through the course of calibration and 
use of a viscometer, the effects tend to be minimized. 
In the calibration of an instrument, heating effects 
result in a negative correction to the calibration 
constant, while in a subsequent determination made 
in the viscometer the effects result in a positive 
correction to the time of flow. Thus, even though 
no corrections are applied, the product of the instru- 
ment constant and the time of flow tends to yield a 
correct result in the viscosity determination. Of 
course this will only be strictly true when the pres- 
sure drop through the capillary is the same for both 
the calibration and the determination, and when the 
liquid whose viscosity is being determined and the 


calibrating liquid have the same heat capacity ar 


temperature coefficient of viscosity. The pertine: 
physical properties of the liquids involved wi 
usually be sufficiently alike, however, to minimiz 
an accumulation of errors in the usual laborato: 
practice of calibrating a viscometer 
graduated as to capillary diameter, by calibrating 
larger capillary with an oil whose viscosity has be 
determined in a smaller one. 


series of 


4.4. Variability of End Corrections With Reynolds 
Number 


The possible existence of curvature in Pt versus 2 
plots for viscometers of the Bingham type, such as 
shown in figures 2 and 3, has not been generally 
recognized. Most observers have confined thei 
calibrations to the higher rates of flow in order to 
obtain more accurate pressure measurements, and 
in addition, there has often existed a disregard ol 
inaccurate estimation of the level head and logarith 
mic head corrections, which may mask the true 
nature of the calibration curve at the lower Reynolds 
numbers. The curvature found in the plots for 
viscometers | and 20, however, is not inconsistent 
with the theory of end corrections and may be 
explained through a consideration of the nature of 
the coefficients m and n in eq (1). 

The value of m is dependent upon the work done 
in accelerating the liquid from essential rest to the 
parabolic distribution of velocities existing after a 
sufficient entrance length has been traversed. This 
work is done not against viscous forces, but only 
imparts kinetic energy to the flowing liquid at the 
expense of the applied pressure. The value of n 
arises in part from the excess work done against 
viscous forces in the entrance length as compared 
with an equivalent length where the distribution of 
velocities is parabolic throughout. At sufficiently 
high Reynolds numbers the kinetic energy of the 
liquid streaming out of the exit end of the capillary 
is dissipated as heat [16] in the enlarged part of the 
viscometer past the capillary. For this condition, 
the phenomena at the exit end will make a negligible 
contribution to the values of m and n. At lower 
Reynolds numbers, however, part or all of the 
kinetic energy may be expended against viscous 
resistance close to the exit of the capillary. The 
value of m for a given Reynolds number will then 
depend upon the difference between the work done 
in acceleration at the entrance end and the amount 
of kinetic energy used in overcoming viscous resist- 
ance at the exitend. The magnitude of n will be the 
sum of the extra work done against viscous resistance 
at the two ends. It follows then that higher values 
of m correspond to lower values of n and conversely. 
Also, when m is constant over a range of Reynolds 
numbers, n should also be constant. Dorsey [5] has 
treated this subject at some length and concludes 
that, for capillaries with square-cut ends, m is zero, 
and » has a constant value of 1.14 up to a Reynolds 
number of 10, while at higher Reynolds numbers, 
m=1 and n has a constant value of 0.57. He 
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wnizes, however, that the values of m and n will 
affected by the geometry of the entrance and 
ends of the bore In the construction of the 
cham viscometers some glass blowing was done 
the ends of the capillaries so that they are not 
cisely square-cut, but evidence slight fire polishing 
the edge of the bore. With this geometry it is 
surprising that these viscometers show an 
ended transition region (approximately R= 10 to 
between the range where m=O and the range 
ove R=60 where m becomes constant. 
For Bingham viscometer 1 the constant value of 
s 1.12, which is in agreement with values reported 
others for capillaries with squarecut ends [3]. 
he lower value of m=0.95 found for Bingham No 
) may be the result of some peculiarity in the shape 
the ends of the capillary bore as a result of glass 
lowing. With the Cannon viscometers M25—1 and 
\[25-2, it was found that m=0 up to at least a 
Reynolds number of 120, which is as high as the 
iibration was carried. This condition is un- 
loubtedly associated with the gradual tapers through 
which the bores of these capillaries expand. It is 
probable that at sufficiently bigher rates of flow 
these calibration curves would also pass through a 
transition region and then assume a constant slope 
corresponding to a value of m which is close to one 


5. Conclusion 


In making the measurements reported here, every 
reasonable precaution was taken so that the results 
would refleet the accuracy of which each viscometric 
method was capable. The viscosity values of the 
four liquids as determined independently by the 
two methods were in better agreement than was 
anticipated, and hence it is believed that both 
methods are relatively free of error. Reviewing 
the techniques and sources of error critically, how- 
ever, it is apparent that the necessary steps in cali- 
brating the Bingham viscometers are more compli- 
cated and contain more possibilities of uncertainty 
than are present in calibrating the Cannon viscom- 
The experimental difficulties in establishing 


eters, 


the calibration curve with water at low Reynolds 
numbers and the uncertain effects of heating in the 
capillary are perhaps the weakest points in the 
method. On the other hand, the relative simplicity 
of the method with the Cannon viscometers makes 
this method singularly attractive and inherently 
more accurate. The greatest difficulty with this 
method seems to lie in adequately correcting the 
head for the difference in the surface tensions of 
water and oil. The procedures described here for 
evaluating this correction are cumbersome and 
time consuming to the extent that less accurate 
estimates of the correction are resorted to in most 
laboratories. It is possible that a redesign of the 
shape of the fiducial bulb of the viscometer would 
simplify the evaluation of this correction. For 
example, a cylindrical bulb with a conical top and 
bottom would present simple geometric shapes for 
which surface tension effects are readily calculable 
by Barr’s method [11]. Any redesign of this nature, 
however, would have to be accompanied by consid- 
erations of the drainage characteristics of the bulb 
as well 
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7. Appendices 


7.1. Calibration of Bingham Viscometer 1 
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7.2. Calibration of Bingham Viscometer 20 with Water at 20° C 
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7.3. Calibration of Bingham Viscometer 5 with Oils at 20° C 
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7.4. Determinations in Bingham Viscometers 
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7.5. Calibration of Viscometers M25 1 and M25. 2 with Water 


st) 


n the times 


sat 40°, 60 


7.7. 


sec 
248. OS 


2A 00 


1110. 
1110. 6 


291. 30 
201. 25 


1085, 11 
1085. 40 


to a hypothetical filling represented by the mean of the 





tnean vl justed MA h SMayih Mach Vavh ! Ah > } 
for ea 
se sec se se 
6. 33 
“3 | 
7 2 Mi. Ht ww 0. OO0OR 0. G0000 0. 00000 0. OOOOH LLL wid 1. woe Ty 1 
42 
182. 72 ‘82.72 in2. 81 oooee ool O0008 ooo! 1. QO00S {A278 1. 3632 “4 
4a. 1 Mal Hi8. 14 +. OOOO ooo + OOO18 ooo Io MA. 16 1. 3620 tw 
" “) » wo wy s. Oooo OOOS6, OO LLL! 1 OOO a ~) 1. a0 l 
4. ¥. | 
1 674.06 O74. ono GO00o ooooo Oeooo mLLLie 4.4 1. 48718 lt 
or | 
4 ] M245 M24 + Oooo? ooo + (OOS +. OOO] 1, OOS 142 48 1 4 : 
419.1 $10.07 110.0 ooo OOOs1 + OOOTB +. OOO 1000 su. 08 1. 4870 1 
5.44 5.44 245.44 aooe2 00056 +. OOODH OOS 1, Oooo 45.44 1.48 mo 
of flow at 20° C for the different fillings are assumed to be largely due to differences in the volume of liquid charged into 
md 80° C Were obtained with different single fillings, the observed times of flow at these temperatures Were adjusted to 
observed times at 20° ¢ 
7.6. Properties of Water 
lem perature i] e "sé 
« Dynesiem cp gem ca 
PU) 72.7 1.0020 0 1. OO88 
“) O. 6541 0. 6482 
“wo wis 74 
x0) 4S wl 
International Critical Tables 
Hardy and Cottington [12], based on » 1.005 ep and converting to 
1.002 ep 
N. E. Dorsey, Properties of ordinary Water substance, table 03, in 
ims per milliliter and converting to grams per cubic centimeter, 
1 1 ml 1.000027 em?* 
Calibration of Viscometers M104 and M105 with Oils 
meat Say/h May/h Saa/h 1+DA/h t ‘ kK Ke, mea 
Se ae ¢ cs 
\ 908 ON +0. 00009 0. 00034 +0). OO144 1.00119 25. 44 2. #42 8. 257210 
j 
& 2M 10 
110 +. OoO08 ooo + oo142 1. OO12 1111.07 1811 & 2M) 
, 201. 28 +. OO008 o0034 ood 1.00119 201.63 2 #42 & 4400 
j 
t noms 
LOS. 26 +. O00 0002 +. OO142 1, 00126 1086.63 | 9 1811 8. 4402 | 
} 


119 


-.R, 


-_- = 


roe SFT tem eee 


SPR PATI 





K 
K 


ih 
Bh 


Heptane M251 
Deo M251 


1) M2 


M251 


M2 
M2 


M25.1 
M25-1 


M 25-2 


M2 


Mio4 
M 104 


Mio 
M105 


7.8. 


nee 
10. 
440. 38 
wm 
Hi. SS 
TSO). 0 
S05. 80 
I" “4 
1655, 28 
AT27. 52 
2 “) 
“le wn" 
Hin. 12 
oT. 22 
v7. 04 
we on 
whl. 1¢ 


WASHINGTON, September 16, 


Results Obtained with Cannon Viscometers 
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Coaxial Radio-Frequency Connectors and 
Their Electrical Quality 


M. C. Selby, E. C. Wolzien, and R. M. Jickling 


The widely accepted manner of evaluating the quality of coaxial radio-frequency 
connectors was in the past limited to a single case, namely, to the condition when the load 
terminating the system was equal to the characteristic impedance of the line on the output 


end of the connector 
the input line 


To broaden this method of evaluation a ‘ 


The quality was expressed as the voltage standing-wave ratio in 


connector” is redefined, and several methods 


are given to find the corrections of various types of connectors with any termination 
Typical results of measurements made by these methods are given for frequencies from 100 


to 900 megacyeles 
quality of transmission lines in general 


1. Introduction 


The accuracy of impedance measurement depends, 
among other things, upon the frequency at which it 
is to be measured. It is generally admitted that 
the frequency range of 30 to about 1,000 Me/s, 
which is under consideration throughout this paper, 
happens to be one of the more difficult ranges from 
the standpoint of impedance measurements. This 
is the transition range where components cease to 
act as lumped-constant elements. Dimensions and 
distributed constants must be considered, particu- 
larly when optimum accuracies are desired. Un- 
fortunately, it is difficult to specify quantitatively 
the accuracy sought because the accuracy of an 
impedance measurement is generally a function of 
both the magnitude and phase angle of the unknown 
One may however indicate its order of magnitude 
as | percent or better for “matched” conditions 
In searching for reliable measurement techniques 
and equipment for this frequency range to meet the 
needs of the National Bureau of Standards, the well- 
known slotted-line method was found the most 
promising After considerable effort, slotted lines 
having the highest available mechanical and elee- 
trical precision were realized. These lines, usable 
down to 30 Me/s, had a uniformity of probe-voltage 
output of 0.25 percent or better over the entire 
working range at all useful frequencies, and a poten- 
tial accuracy considerably higher than any other 
available equipment. To make use of that potential 
accuracy, it was necessary to minimize the errors 
introduced by various connectors required between 
inknown impedances and the uniform section of the 
precision slotted line. These errors, usually negli- 
ible with equipment of lower accuracy, could no 
longer be disregarded for full utilization of the 
quality of the new lines. This led to an investigation 
of methods of determining connector quality in 
reneral with results presented below. 

The determination of errors introduced by con- 
necting elements has been a problem for a number of 
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Application of these methods are also indicated in determining the 


vears, particularly when used with measuring or 
controlling equipment, e. g., attenuators, terminating 
loads, matching indicators, power monitors, bridges, 
etc. In most cases these connecting elements con- 
sist of rigid coaxial structures, referred to as “con- 
nectors,”’ with or without a short solid-dielectric 
cable or air-dielectric rigid transmission line. The 
characteristic impedances of the components of 
these connecting elements usually differ from each 
other by several percent. In addition, the connectors 
may have a nonuniform internal structure and dis- 
continuities of their own. The losses are usually 
negligible and are assumed to be so throughout this 
treatment. 


2. Definition of Connector and of Its Quality 


The tendency has been to interpret a “connector” 
as the physical body of a connecting element, such 
as an adapter, a balun, a taper, a jack, or a plug, by 
itself or in combinations without regard to discon- 
tinuities located immediately at the terminating 
planes. Thus, some indefinitely located plane be 
tween a type N plug and jack was taken as a reference 
plane for all purposes, e. g., the plane where a short 
circuit was to be located, the plane where the un- 
known impedance, voltage, or current was to be 
measured, ete. When accurate knowledge of im 
pedance is desired, this interpretation of a connector 
is not always acceptable, because for TEM propaga 
tion (considered here) impedance Is basically defined 
and treated as the ratio of voltage to current or of 
electric to magnetic fields in a cross sectional plane 
of a uniform distributed-constant system. There 
are no known commercial connectors meeting the 
requirements of such a system: all have appreciable 
discontinuities affecting the field distribution in a 
manner difficult to compute, predict, or measure. 

The term “connector” as used here includes all 
lossless circuit arrangements used to interconnect 
two transmission lines of any configuration and 
characteristic }impedance. An example of such a 
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is shown in figure 1. The transmission 
assumed to have uniform distributed 
constants up to the connector To complete the 
definition of a connector, it is necessary to fix the 
minimum allowable distance between its terminating 
planes and its adjacent discontinuities (e. g., h, fig. 2). 
This distance depends on the particular configuration 
of the transmission lines and the attenuation rate of 
the higher modes generated at the discontinuities. 
For coaxial lines, A (fig. 2) should be equal to or 
larger than the diameter of the outer conductor [1}.' 
When used in this sense, connectors can be readily 
treated as four-terminal networks and their constants 
defined and measured. 

Another accepted convention is to indicate the 
quality of a connector by the voltage standing-wave 


connector 


lines are 


' Figures in brackets indicate the literature references at the end of this paper 





ratio (VSWR) introduced on a uniform transmissi: 
line feeding the connector when the latter is tern 
nated in a load equal to the characteristic impedan 
of the line. This convention can only be used 
estimate the maximum error in a measured VSW) 
caused by connector reflections. It does not furnis 
necessary corrections fer all impedance valu 
The approach in this paper is to treat the connect: 
as a component part of a transformer; the figuri 
of merit of this transformer are two transformatio 
corrections for the components of its terminatir 
impedances. The problem is to establish relatio: 
ships between the connector input and output 
impedances or admittances. Admittance relatio: 
ships prove to be more convenient for computatio: 
simplicity. 

Let, in figure 2, a uniform transmission lin 
having a characteristic impedance Z , be joined 
with another uniform line, Zm, by means of a con 
nector terminated in the planes A’, at its input end 
and Sp, at its output end. Also, let 7) be the location 
of a voltage node when a short circuit is placed at S 
It can be shown [2] that simple expressions inte 
relate admittances on opposite side of the con 
nector, e. g., the admittances Y,=—G,+)B, and 
VY. =Gy + 7By at planes A and A’, where A and 
A’ are effectively an integral number of half wave 
lengths apart. Two planes in a transmission-line 
system are effectively an integral number of half wave 
lengths apart when a voltage node placed in one of 
these planes (closer to the load) will cause a voltage 
node to appear in the other plane irrespective of 
what the physical distance is between them. A sup- 
plementary line section, /,, is added to fulfill the 
above requirement. The relationships of these 
admittance components (fig. 2), as shown by Oliver 
[2], are given by 

Gy ara, (1 
and 
Bi =a?B,+a,b;, (2 
where a,’ is a dimensionless constant, and a,b,’ is an 
additive constant in mhos. Similarly, 


a", a 2G, (3 
and 

Bi, =a By +-ayb,. (4 
If desired, one can treat the connector shown in 
figure 2 between A’ and S, (minus the supplementary 
line section /,) as a transformer by expressing the 
input admittance in terms of the output admittance 

From well-known interrelations 


oa g,(1+ cot? Bol) (5) 
o'* (cot Bal,—b,)?+92 


_b, cot? Bol,+(1 —gi— bf) cot B.l.—b, 
(cot Bol,—b,)? +95 


(fh 


b, 


where: 8, is the phase constant and the admittance 
components are normalized with respect to Yip. 
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Che termination admittance of the connector may 
, be obtained, using various charts and diagrams 
thus avoiding computation of eq (5) and (6) 
wever, there is no great advantage in knowing 
admittance at the plane S) in preference to that 
any other plane. The additional computation 
k involved is therefore seldom justified. 
lhe definition of a “connecting transformer’ may 
be introduced as follows: 


\ connecting transformer is any linear bilateral 
passive four-terminal lossless network terminated at 
both ends in uniform line sections; the latter are 
electrically and mechanically identical to the corres- 
ponding adjacent lines of the system. When a short 
circuit is placed at either end of the connecting trans- 
former, a voltage node appears at its other end 

ius the connecting transformer consists of the 

nnector and 1 or 2 supplementary line sections. 
llowever, there is generally no practical advantage 
of having more than one supplementary section; the 
lisadvantage is, of course, the need of two (instead 
f one) admittance transformations when one is 

iterested in the actual input and output admittance 

{f the connector. A connecting transformer, as 
defined above, therefore consists of the connector 
proper and one supplementary line section at its 
output or at its input end. In figure 2, 2, is the 
input-end supplementary line section of the con- 
necting transformer correlating Y; and Y, through 
a set of corrections aj and a,6,. If the output 
supplementary line section /, is chosen, a set of cor- 
rections, @ and a,6,, is required, and the correla- 
tion between ¥", and Y, is established. It is con- 
venieat to use /, in preference to /, when the quality 
of connectors is studied. On the other hand, /; is 
used when the connector is an integral part of an 
admittance measuring system used to determine Y, 
with optimum accuracy. The justification for these 
preferences will be further clarified below when the 
methods of obtaining the transformation corrections 
are discussed. Special problems may indicate a 
different choice of supplementary sections and may 
justify the use of two sections simultaneously (i. e., 
one at the input, the other at the output). How- 
ever, the two major groups of problems anticipated 
at present are impedance and connector measure- 
ments; therefore the discussion below will be limited 
to connecting transformers using either /, or /,. 

The interpretation of a’ and ad’ in the case of 
impedance measurements is straightforward. For 
connector measurements the interpretation is as 
follows: With a perfect connector (referring to fig. 2) 
normalized admittance ¥",/ Yq; will always be equal 
to the normalized Y,/Yq; otherwise the values of 
a and ab’ will indicate the errors introduced by the 
connector discontinuities, i. e., the relative quality of 
various connectors for all values of admittances. 
For a connector mounted on the panel of an enclosed 
network like an attenuator, voltage generator, etc., 


one may simply consider the connector as a part of 


the enclosed network with the reference plane at 
\'. Otherwise, 5 ger reference impedance-planes 
ke Sy and A may be assumed inside the network. 
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To summarize, the definitions of connector, con- 
necting transformer, and connector quality, as 
applied here, are as follows: (a) A connector is a 
lossless linear bilateral passive four-terminal net- 
work terminated at each end in a uniform line section 
of a minimum length required to attenuate higher 
modes likely to appear at these ends. These line 
sections are identical electrically and mechanically 
to the respective lines of the interconnected system 
(b) A connecting transformer is a connector as 
defined in (a), terminated at one end in a supple- 
mentary lossless line section. This supplementary 
line section is less than half a wavelength long and is 
identical electrically and mechanically to that at the 
connector end. When a short circuit is placed at 
either end of a connecting transformer (with r-f 
power fed from the opposite end), a voltage node 
will appear at the other end of the transformer. 
With perfect connectors a connecting transformer 
is an ideal transformer; the normalized input admit- 
tance components of this transformer are equal to 
its normalized output admittance components. 
(Normalization is with reference to input and output 
lines, respectively.) (c¢) With imperfect connectors 
one must apply corrections to obtain accurate 
termination admittances. These ‘‘transformation 
corrections” a? and ab’ indicate the quality of the 
connectors. The simplest form of the corrections 
are given in eq (1) to (4). These transformation 
corrections hold only for the chosen combination of 
connector and supplementary line section. 


3. Determination of Connector Quality 
3.1. Analytical Background 


The inadequacy of specifying the VSWR of a 
“matched”? connector as its quality index has been 
pointed out above. For systems of equal char- 
acteristic impedances the effective characteristic 
impedance of the connector is frequently taken as the 
square root of the product of the measured open and 
shorted input impedances (i. ¢., Zo=y ZocZse). This 
vields, strictly speaking, an imaccurate parameter, 
because the approach assumes infinite and zero 
impedances at open- and short-circuited conditions, 
respectively, and uniformly distributed constants, 
none of which are true in practice. In addition, this 
approach neglects likely discontinuities at the 
termination planes of the connector as the latter is 
conventionally referred to. Therefore, the value 
of Z so obtamed cannot be used as a simple trans- 
formation factor. A reliable way to find quantitative 
values of the transformation corrections of connectors 
is the one employing the node-shift method [2, 3, 4, 5}. 
Frequency-variation methods can be used only to 
find qualitative data. 

The node-shift (n-s) method is well described in 
the references cited and consists briefly in measuring 
corresponding shifts of voitage-minima at both ends 
of the connecting transformer. The original appli- 
cation of this method involved a less efficient and 
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less accurate handling of measurement data, re- 
quiring the determination of tangents to curves of 
uncertain regularity. An improvement, introduced 
by Oliver, largely eliminated this difficulty, because 
the transformation corrections are obtained from a 
curve that is analytically shown to be a straight 
line; measurement errors can be readily detected 
in the plot of a straight line. Techniques developed 
for the practical application of this method in systems 
employing equal or unequal characteristic impedances 
are described in section 3.3 Several variations in 
these techniques may be used, depending on condi- 
tions and requirements 

The essential analytical underlying this 
method of approach are as follows: Let a, 6, ¢, d 
represent, at one frequency, linear parameters [2] 
of a four-terminal network equivalent to either one 
of the two connecting transformers shown in figure 
2. As the elements of the network are assumed 
passive and linear, and a zero termination impedance 
is transformed into a zero input imedpance, 


steps 


ad be l (@) 
e=@0 (8) 

ye Z, 
< 3 9 
4\= At abZ, ) 
Yi=a*y\+ab (10) 


In the last two expressions a and 6 are the parameters 
of the transformer between ) and Y, (i. e., between 
planes 7) and S,) 

\s pointed out above, a voltage node at S, will 
cause a voltage node at 7). The plane 7, may be 
located at the input to the connecting transformer or 
at any other plane an integral number of half wave- 
lengths away from it. Let the voltage node now be 
moved a distance s from Sy) by means of a movable 
short circuit. This is equivalent to varying the re- 
actance of Z, by a given amount. The node at 7, 
will therefore shift in turn by some distance ¢, not 
necessarily equal to s, as shown in figure 3. Here the 
transformer input plane is shown at A and 7, is 
chosen a half wavelength away from A; this is neces- 
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sary to permit measurements along a uniform | 
With the nodes in these positions, 


section. 
Z,’ Zo tan §,f 
Z, == jZe tan Bos 


where 8 
lines. 


2nd is the phase constant of the respect 


Using eq (9), 


ya tan Box 
7 tan £,l ; ; ’ 
Loa a*- JabLZ o tan Bs 


since the network is lossless, 6 is an imaginary quay 


tity [2] and one may substitute —6’ for 7b. Rew: 
ing the terms, 
sk , 
cot B,t=+" (a? cot Bos—ab’Zo), 14 


therefore, 


cot B,t—cot Bos (> a*—1 ) cot Boxs—ab’Zo,. (15) 


One can thus obtain a number of cot 8, values 
corresponding to chosen cot #« values and the plot 
of (cot Bt—cot Bs) versus cot Bs will be a straight 
line with |(Z»,/Z 2)a’— 1] as its slope and (—ab’Z 
as its intercept on the (cot 8,f—cot Bes) axis. As 
Z, and Ze are assumed to be known, a’ and ab’ can 
be determined for use in eq (1) to (4). In addition, 
relationships for impedance components can be 
derived to show how the load impedance, R+-7.X, of 
a slotted line, for example, is determined from its 
measured value, R’+j7X’, when the transformation 
corrections of its connecting transformer are known 
Thus, from eq (9), again writing — 6’ = jb, 


a(R’ 4 px’ ) 


BT IN TR jab (R'+ 5X) (16 
or 
R — oR = bes ‘ 
1+2ab’X (ab) (R724 X") 
, a*LX’ L@bh’(R’?4 X”) 
XA 1+ 2abh’X’ j (ab’y(R” j xy") 18) 


It may be of interest to note at this point the inter 
relation between the parameters of the same connec! 
ing transformer in the forward and reversed dire« 
tions. A case may arise where it is necessary to 
reverse the transformer end for end. For exampl 
because of external structural irregularities, one ma) 
wish to connect the transformer into a circuit in a 
direction opposite to that used for the measurement o! 
its constants. The termination impedances mus! 
then be corrected by a new set of parameters, because 
the system is not necessarily symmetrical. There is 
a simple correlation between the “forward” an 
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ersed”’ sets of corrections, obviating the necessity 


dditional measurements. It can be shown [6] 
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re the subscripts f and r refer to the forward 
reversed conditions, respectively 

et the intercept of the straight-line plot measured 
he forward direction be n; then 
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also the slope obtained in the forward direction 


then 
a?==— (m+ 1 22) 
TY reversed svstem 
P bh’ a,b’ n 
a,b 23 
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Kquations (22) and (23) may be used to compute 


the output impedance of a connecting transformer 
The impedance at the input end of the connecting 
ransformer is assumed to be known. 


3.2. Connectors in Tandem 


There is frequently the need of measuring the 
quality of two connectors in tandem im order to 
‘btain the quality of one of them. A_ typical 
example is as follows: Precision slotted lines of 
special mechanical and electrical design are fre- 
queatly constructed at a considerable cost; special 
liameters are chosen to obtain optimum quality of 
he slotted line. It is therefore necessary to employ 
lapters, tapers, etc. to connect this line to conven- 
onal connectors of various diameters. There is a 
simple relationship between the transformation 
orrections of a connecting transformer composed of 
vo transformers in tandem and the corrections of 
he individual transformers 

In figure 4 leta precision slotted line be connected 
hrough connector | and connector 2 to an unknown 
ad in plane Sy at the ead of connecting transformer 
2. The supplementary line section of transformer | is 

a section of the precision slotted line. From the 

1) to (4), 

a,b,, (24) 


dy2b4,—a2a,b; 


here subscripts 1, 2, and 12 correspond, respectively, 
connecting transformers |, 2, and 1-2, and where 


mes Sy, and 7) coincide. 
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It is thus seen that one can obtain the individual 
transformation corrections of transformer 2 from 
the corrections of the combined transformer 1-2 
This is a great advantage, as it allows the use of 
slotted lines of the highest precision to accurately 
measure connector quality and unknown impedance, 
even when the latter do not match the line mechan 
ically or electrically 

The significance of the magnitudes of the trans 
formation corrections a? and ab’ may be demonstrated 
as follows: Equation (1) indicates that the error in a 
measured conductance is a function of a*, whereas eq 
(2) shows that the error in a measured susceptance 
is a function of a, 6’, and of the magnitude of that 
susceptance. From eq (2) the fractional error ta 
this susceptance is 


BB ' ab’ ~— 
= a”) (20) 
B B 
or 
r—z’ ' 
, ab’r+(a 1), (26) 
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where B l/r, B’ 1/x’; x and 2’ are the respec- 
tive equivalent shunt reactances of the unknown and 
measured impedances at the ends of the connecting 
transformer. 

Equation (26) thus gives the errors in the measured 
shunt reactances, which can be plotted versus 2 as 
straight lines with shunt capacitive reactances 
equivalent to (—ab’) mhos as slopes, and (a4?—1) as 
intercepts. Representative values usually encoun- 
tered in practice at frequencies to 1,000 Me/s are 
shown in figure 5. 


3.3. Techniques and Procedures 


The node-shift method of connector-quality deter 
mination may utilize one of the following reactors, 
depending on the problem and facilities at hand: (a) 
An adjustable stub, (b) fixed line sections, (¢) a 
slotted line, or (d) a perforated line. Steps of pro- 
cedure and precautions that are common to the use 
of all of these reactors will first be discussed; then 
features of each device and its associated technique 
will be described individually. 


a. Equipment 


In general, the equipment requirements are these: 
The known reactor, connected as shown in figure 3 
at Sy (the end of the unknown connecting trans- 
former) must be lossless, uniform, and adjustable or 
variable over at least half a wavelength. The junc- 
tion at S, must be free of discontinuities; thus, for 
best results, the cross-sectional dimensions and di- 
electric material of the reactor and the connecting 
transformer must be the same at this plane. In 
addition, this uniformity must be maintained on the 
input side of So for an axial distance equal to at least 
the diameter of the outer conductor. Sometimes, 
as mentioned in section g, a small amount of mis- 
match may be tolerated without invalidating the 
measured data. Good mechanical alinement and 
contact at this plane are essential. All discontinui- 
ties should lie between Sy) and A, the latter being a 
plane of reference on the measuring line, effectively 
an integral number of half wavelengths from Sp. 


b. General Procedure 


The measurement procedure consists in first plac- 
ing a short circuit at S) and locating 7), a position on 
the measuring line where a voltage node occurs, one 
or more half wavelengths from AK. Whenever pos- 
sible, the first 7, beyond K should be used so as not 
to introduce extra losses into the measuring system. 
As the short circuit is moved from S, to S, a distance 
s, the voltage node on the measuring line moves from 
T, to 7, a distance ¢. From 5 to 10 corresponding 
positions of S and 7 are recorded. This number of 
points has been found desirable to identify any ir- 
regularities in the data. Values of s are chosen so 





that the distribution of cot B.s will be fairly even 


spaced from —2 to +2. These limits are select 
because usable data, for the rapidly changing cota: 
gent function outside this interval, require equipmer 
precision which is usually beyond practical achiev: 
ment. When the magnitude of the s values is lin 
ited by the usable length of the reactor, lengths of 
corresponding to cot 8,8 values from 0 to +2 may | 
adequate. Each of the S and 7 nodes are located bh 
determining the midpoint between two positions o 
equal detector response one on each side of the mini- 
mum. However, when structural defects in th 
slotted line produce nonuniformity in the voltage 
distribution along the line, errors may be caused i 
the determination of the nodal positions. To mini 
mize these errors, it has been shown [7] that the opti 
mum equal-response voltage in the determination of 
a nodal position should be approximately 3 db 
(within +0.5 db) above the minimum voltage when 
the VSWR is above 7. In addition to structural 
defects, the presence of losses in the measuring SVs- 
tem may introduce errors in making accurate nodal 
measurements. When such is the case, one may 
minimize these errors by following the procedure 
outlined later in section g. 

The difference, cot 8,f—cot Bs, is then plotted 
versus cot 8s and the resulting straight line inter- 
preted as shown in section 3.1. The values of @ 
and ad’ thus obtained are the transformation correc- 
tions of the unknown connecting transformer at the 
test frequency. 

As the magnitudes of a? and ab’ are functions of 
the relative positions of the discontinuities in the 
connecting transformer, and as these positions are 
scattered and unknown, the effect of these discon- 
tinuities may cancel at some frequencies and increase 
at others. Therefore, the number of frequencies at 
which measurements must be made in order to find 
the quality of a connector over a given frequency 
range will be determined by the extent of the fre- 
queney range and the accuracy desired. Measure- 
ments indicated that the transformation corrections 
for commercial connectors varied rather slowly with 
frequency up to 1,000 Me/s, and that measurements 
at 100-Me/s intervals in this range would be sufficient 
to determine the frequency characteristics of a given 
connector. 
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necting transformers in tandem. 
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c. Connectors in Tandem 


When the measuring line does not match the un- 
wn connector in size or characteristic impedance, 
both, an adapter is required between them. Then 
s necessary to separate the corrections of the 
cnown connector from those of the adapter ac- 
ding to the procedure described above in section 
In this situation Sw» (see fig. 6) must be so 
ed that a voltage node appears at Sy. This 
ulition can be met by first placing a short circuit 
Sy, the output end of connecting transformer No. 
and by locating a corresponding 7), on the meas- 
ng line. Then the short circuit is replaced by the 
known connector and adjustable reactor, and the 
tter is varied until a voltage node again appears at 
Thus a voltage node is also assured at Sy. The 
lues of s and ¢ are then obtained as described 
viously. In some cases the voltage node at Sy 
inot be moved continuously because of the type 
variable reactor being employed or because con- 
cting transformer No. 2 has a fixed length. The 
rocedure to place a voltage node at Sy, under such 
maditions will be deseribed in section h. 


d. Curvature in Plot of Data 


It was occasionally observed that the plot of meas- 
wement data substituted in eq. (15) resulted in a 
ive rather than a straight line. The following 
precautions were found to prevent such curvature. 

The phase constants of the measuring line and the 

reactor must be known to an accuracy of +0.01per- 
cent. This required that the corresponding fre- 
quency stability be maintained to an accuracy of at 
least +0.01 percent. (2) The positions of S, and 7, 
should be measured, for example, to +0.1 mm at 300 
\e/s for lines with dielectrics having constants of 1 
to3. The presence of any appreciable discontinuity 
at these planes, such as a poor contact at Sy, may 
cause an excessive error. (3) When the voltage 
distmbution near a node is nonsymmetrical because 
of losses in the measuring system, proper precautions 
must be taken to locate the positions of S and T cor- 
rectly. These precautions are described in section g. 
1) The accuracy of finding voltage-node positions 
may be improved by keeping the input voltage high 
and the probe penetration deep (for low input volt- 
ages), probe reflections being negligible near voltage 
nodes, by filtering out harmonics from the input 
voltage, and by operating probe output indicators 
at the most sensitive positions of the meters and cir- 
cult action. Means of minimizing errors caused by 
structural defects in slotted lines were indicated 
under section b. 


e. Use of an Adjustable Stub as a Reactor 


This type of reactor is the most obvious one to 
select for use in n-s measurements of air-dielectric 
elements and is quite simple to operate. A 5-ft 50- 
ohm stub of this type was found particularly useful 
in measuring %-in. connectors in the VHF range. 


To determine values of s, relative stub positions 
may be measured with the required precision by 
means of large vernier calipers, gage blocks, or a 
carefully engraved scale attached to the adjustable 
stub. An advantage of this type of reactor is that 
one stub of sufficient length can be used to obtain 
convenient values of s over a wide frequency range 

In the construction of the stub it is essential to 
use the least amount of solid support for the center 
conductor near the input end. Whatever support 
may be required should be small in mass and have 
the lowest possible dielectric constant. When it is 
possible to support the center conductor by the con- 
nector being measured, the omission of any additional 
support is desirable. Another critical construction 
feature is that the movable shorting device should 
locate a nearly lossless short cireuit in positively 
defined planes. For small-sized connectors it is very 
difficult to build a stub long enough for use in low- 
frequency n-s measurements without an objection- 
able sag in the center conductor. An adjustable 
stub is unsuitable for measurements of solid-dielectric 
connectors, because of unavoidable discontinuities 
of the junction of the two dielectric media 


f. Use of Fixed-Line Sections as Reactors 


The use of accurately constructed line sections for 
reactors simplifies ns measurements. First, a short 
circuit is placed at S); this is then replaced by 
shorted fixed-line sections of different lengths corre- 
sponding to selected s values Alternatively, sec- 
tions of equal length open at both ends may be con- 
nected in tandem; a short circuit is then placed at 
the output end of the combination. In both cases 
it is essential that a good mechanical contact exist 
at S, and that no discontinuities be introduced in 
this plane In using combined sections, these condi- 
tions must also be maintained at the junctions. 
Otherwise, a curvature or excessive random scatter 
will result in the final plot of the data 

Although the procedure followed with this type of 
reactor is simple, a number of disadvantages are 
involved. First, the careful matching and alinement 
referred to above requires precision machining. 
Thus when n-s measurements are to be made over 
a wide frequency range, numerous costly precision 
sections are required. Second, if an adapter Is 
required between the unknown connector and the 
measuring line, it may be awkward to use fixed 
sections as reactors (particularly at low frequencies) 
because of the requirement that the ends of the 
unknown connecting transformer be an effective half 
wavelength apart (see section 3.2). Third, although 
reactors of this type can be constructed to match 
any present-day air-dielectric connector, this can not 
easily be done for solid-dielectric components because 
it is difficult to determine with sufficient accuracy 
the phase constant of solid-dielectric sections In 
coaxial lines having solid insulation the phase con- 
stant is very closely equal to 2f(k)'’?/(2.998 & 10°) 
radians/meter, where & is the dielectric constant of 
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the insulation. 
be found in the dielectric constants of commonly 
used insulators, depending on the batch or manufac- 


Variations from 1 to 3 percent may 


turer. In addition, these constants may be affected 
by a change in volume (i. e., density) caused during 
fabrication steps and by the application of excessive 
heat. The phase-constant accuracy is also limited 
by the accuracy (+ 0.1 percent) to which the capacity 
of line sections can be measured, by effects of fre 
quency on the dielectric constant, and by dimensional 
tolerances of the conductors. Thus it appears that 
the best prospect of obtaining the required accuracy 
of +0.01 percent in the phase constant is to measure 
the electrical lengths directly in the finished line 
sections. For example, a section of line can be 
resonated by varying the frequency within the range 
in question and observing the output of a fixed probe 
located at a precisely known distance approximately 
a half wavelength from the shorted end. Assuming 
that the dielectric and dimensional uniformity was 
maintained, this line section can then be cut into a 
number of short sections. This fixed probe method 
suggests the application of a perforated solid -dielec- 
tric line; it also suggests an accurate way to measure 
the dielectric constant of a coaxially fabricated 
material, 

A type of reactor combining the features of the 
adjustable shorting stub and _ fixed-line 
might be used to alleviate some of the disadvantages 
mentioned above. First, an air-dlielectric stub is 
required, in which the shorting device can be removed 
from the coaxial line. Coaxial sections of any de- 
sired solid insulator are then inserted in the coaxial 
line, and the shorting device is replaced so as to 
produce a short circuit at the end of the newly 
formed solid-dielectric coaxial line. It is necessary 
that each of the insulator sections be accurately 
machined to fit the coaxial line and to have end 
planes perpendicular to its axis. The lengths of 
these sections depend on the desired values of cot 


sections 


8.x. Advantages of this type of reactor are (a) 
the conductors are not broken: (b) numerous s 


values are readily available, because it is simpler 
to machine any desired number of insulators rather 
than to construct an equal number of individual 
line sections; (c) the terminiating plane of the 
unknown connecting transformer  (solid-dielectric 
type) can easily be placed at any approximate point 
in the reactor, when required; (d) any number of 
different dielectric materials may be used with the 
same air stub for various desired Z» values: (e) the 
phase constants may be determined with an auxiliary 
fixed probe, as described above. 


g. Use of a Slotted Line as a Variable Reactor 


In employing a slotted line as a known reactor, an 
auxiliary line stretcher or adjustable stub is connected 
at the output end of the slotted line to move a voltage 
minimum along the slotted line where it is detected 
by a probe (see fig. 7). The method recommended 
for accurately locating this voltage node is the same 
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as that described in section b. It is not essentia 
that the voltage node be moved to the exact S 
location previously selected; it is possible to locat 
the position of the voltage minimum very accurately 
and thus know the corresponding value of s 

At higher frequencies the presence of losses in 
the measuring system may make accurate measury 
ments of S and 7 more difficult. As the voltag 
curve in the region of a minimum closely approaches 
a parabola, the voltage nodes may be located ac 
curately by joining the midpoints of three horizionta| 
chords through each selected trough of the voltage 
standing-wave [8]. However, this procedure is 
quite tedious; in most cases the following simple: 
method may be employed, with results of sufficient 
accuracy; (a) place a short circuit at Sy and find the 
true 7); (b) connect the slotted-line reactor at S 
this will introduce losses (especially when a_solid- 
dielectric line is used), reducing the VSWR and 
upsetting the symmetry of the voltage distribution 
in the measuring line particularly; (c) locate a nod 
at 7, as follows: Adjust the auxiliary variable reactor 
to obtain a node near 7), using a voltage increment, 
AE,, about 3 db above the minimum; readjust the 
variable reactor to move the node toward the desired 
location; repeat these steps until a node is obtained 
exactly at the true 7); (d) then locate the correspond- 
ing minimum for S;, (an effective half wavelength 
from S, toward the short circuit). Use a certain 
AE, approximately 3 db above the voltage minimum 
(e) Finally, locate the selected S and corresponding 
T positions, using the same AF, and AF, in each case, 
as was used to locate the Sj and 7), respectively 

Relative merits of this technique are: (a) This 
type of reactor is useful wherever a slotted line is 
available to match the connector being measured 
The technique was found particularly applicable in 
determining the quality of \%-in. connectors, the size 
in which solid-dielectric slotted lines of fairly good 
precision have been made. (b) Only one slotted line 
is needed as a reactor for n-s measurements over a 
wide frequency range, e. g., from 100 to 1,000 
Me/s. (c) Because of the precision with which 
voltage nodes can be located, the velocity of propaga- 
tion, and hence the phase constant of this type of 
reactor, can be very accurately determined. (d) The 
use of this reactor is more convenient than the others 
for the n-s measurement of connectors in tandem, 
because of the ease with which S» can be located 
whenever it must be placed in the reactor to meet 
the half-wavelength requirement referred to 
previously. 
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turn, this technique possesses a few disad- 
Measuring the s values in a slotted line 
omewhat more time consuming that adjusting 
n. as described in the previous methods A 
ted line is considerably more expensive than 
ctors of the other types for use over the same 
jueney range. Because of the slot, it is impossible 
natch perfectly a slotted line to a closed line, a 
promise necessarily being required either in the 
duetor dimensions or in the characteristic im- 
ances However, the use of a slotted line as a 
ctor is practical if the resulting discontinuity at 
junction of the line and the unknown connector 
kept as small as possible. For example, satis- 
tory data were obtained by using such a slotted 
to measure the quality of adapters on precision 
tted lines 


tages 


Use of a Perforated Line as a Variable Reactor 


Although applicable to making n-s measurements 
connectors in all sizes, this technique is especially 
commended for use in determining the quality of 
rve-sized connectors because of the economy made 
possible in the construction of the variable reactor, 
\s an example, a 50-ohm 3\-in. perforated air- 
iwlectric line, shown in figure 8, was constructed for 
the quality measurements of connectors and lines 
having that diameter and characteristic impedance. 
The line consisted of a 4-ft brass outer conductor 
vith a duralumin center conductor supported at the 
output end by a Teflon dise \ in. thick. At the 
nput end no dise was used; the center conductor was 
supported by the connector being measured so as to 
minimize the discontinuity at the junction of the 
connector and reactor. Stock tubing was used for 
the conductors. Probe holes, \ in. in diameter, 
spaced 1 in apart, were bored in the outer conductor. 
This spacing permitted the choice of convenient s 
values over a wide frequency range. Surrounding 


each hole was a milled “‘flat”’ to help center the probe. 
The probe was a stiff copper wire, 0.020 in. in diam- 
eter, soldered to a BNC connector, and was shielded 
by a short cylinder that fitted the probe hole. An 
external adjustable stub was used to tune the probe 
see fig. 9). As with the slotted line, a _ short- 





Fiaure 9. Perforated line 
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Figure 8&8 


Details of a perforated line used as a variable 


reactor in node-shift measurements 


circuited line stretcher was used at the output end to 
move the voltage node along the perforated line To 
facilitate placing the voltage node at a preselected 
hole, a line stretcher that can be moved by very 
smali increments ina positive fashion is desirable. If 
difficulty is experienced in this operation (as indi- 
cated by a lack of reproducibility), the following 
method may be used: (a) Shift the voltage node 
slightly to one side of the selected hole Let the 
probe output voltage at this hole be V, and the 
position of a corresponding voltage node on the 
measuring line be T, (b) Move the voltage node to 
the other side of the selected hole so that the new 
probe voltage V, equals V\. Let the position of the 
corresponding voltage node on the measuring line be 
T,. (c) Then the true 7 for a voltage node placed 
at the desired hole will be midway between 7, and 
T,. This procedure assumes that the voltage curves 
at S and 7 are not distorted by excessive losses 
This is not always true at 7. If such is the case, the 
precautions discussed in section g must be observed 
in finding 7. In all cases, when the movable short 
circuit is being adjusted, it is advisable to have the 
source impedance equal to Z, of the line connected 


used as a variable reactor to measure the corrections of a transformer. 






to it; this assures a constancy of voltage amplitudes 
along the system. 

When it is necessary to use an adapter between the 
unknown connector and the measuring line, the 
corrections of the two connectors must be separated, 
as previously explained. Because of the fixed probe 
positions of the perforated line, a special step must 
usually be employed here to determine the exact 
frequency at which the unknown connecting trans- 
former is effectively equal to a half wavelength. A 
short circuit is first placed at the end of the adapter 
to the measuring line, and three values of 7), are 
found for three different frequencies around the 
estimated frequency at which the unknown connect- 
ing transformer meets the length requirement. 
Then the short circuit is removed, and the unknown 
connector and perforated line are attached, the probe 
of the perforated line is placed in a selected hole, and 
three values at 7) are found for three similar fre- 
quencies With these data two curves of frequency 
versus 7) are plotted. The intersection of these 
curves determines the desired frequency. If the 
effective discontinuity of the adapter is small (so 
that the transformation corrections remain 
tially the same for small displacements of S),), it is 
permissible to use an approximate frequency, e. g., 
the integral number of megacycles per second 
nearest the intersection of the above curves. Then 
the original corrections of the adapter, predetermined 
at aslightly different frequency, can still be used in the 
separation procedure. When it is necessary to know 
the desired frequency more precisely, the above pro- 
cedure is repeated over a more narrow frequency 
range bracketing the frequency located above. The 
frequency used in this technique must be known and 
kept stable to 0.01 percent. The corrections of the 
measuring-line adapter may, of course, be remeasured 
at the new exact frequency if optimum accuracy is 
to be assured 
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TABLE 1 


Connecting transformer 


Number Components 


e-in.-diam line; taper to 2! 


line 


e-in.-diam line; taper to 15s-in.-diam solid-dielectric slotted line 


Adapter UG-57 

Adapter UG-5 

Adapter UG 

Adapter UG— 

Adapter UG-1 

Adapter UG- 

Adapter Two GR-&S74-QNP * 
Connectors UG-21B/U, UG-23B/T 
UG-21B/U, Dem, RG-WU, UG-21B/t 


Step adapter type N to 74-in. diam.* 
7s-in.-~<diam, 140-em partial-air-dielectric line 
?s-in.-diam, 150-cm Teflon-bead line 


34e-in.-diam, 140-cm partial-air-dielectric line 


3\s-in.<liam, 1%)-cm air-dielectric line 


, UG-23B/U mating connectors were used; 4=1.5 in. 
mating connector at input end 


¢-in.-diam solid-dielectric precision slotted 


4. Typical Measurement Results 


In all the techniques described above, a precisi 
slotted line has been used to measure the inp 
reactance variations of the connecting transfora 
corresponding to known output-reactance variatioy 
Obviously, impedance-measuring devices other th 
slotted lines may be used for this purpose; e 
bridges, Q-meters, instruments emploving direction 
couplers, admittance comparators, ete. lowev: 
present experience indicates a need of precision 
measuring input-reactance variations considerab| 
greater than that provided by currently availab! 
commercial instruments other than slotted lines 
Also, when these other impedance-measuring devices 
are employed for this purpose, the steps of proc 
dure in the n-s method must be modified. A dis 
cussion of the use of such devices is reserved for 
another treatment of the subject. 

The results of some n-s measurements of typical 
connecting transformers, using a precision slotted 
line as the measuring instrument, are compiled in 
table 1. The transformation corrections given do 
not necessarily represent average values for the 
particular types of commercial units listed, but 
rather demonstrate the typical magnitudes of these 
corrections. Each connecting transformer, of course, 
terminates in uniform line sections and, when re- 
quired, includes mating connectors. For type N 
connectors and adapters the mating connectors are 
assembled from UG-—21B/U plugs and UG—23B/U 
jacks, combined with 1%-in. uniform sections of 50- 
ohm Teflon-dielectric rigid line. For a particular 
connector the transformation corrections vary, de- 
pending on the position of the connector along the 
transmission-line type of transformer. The 
rections for the cable connectors and adapters, Nos 
3 to 11, in table 1 represent the parameters of each 
transformer, with the test unit placed at the input 


cor- 


Transformation corrections of typi al connecting transformers 


Frequency Variable reactor 


Micromhos 
410 
20 


—f) 
430) 
360 


160 {6-in.-diam solid-dielectric slotted line 
110 7o= 49.9 ohms 
yw) 
—130 
0 
~ 250 
—2 
—50 
680 





— 90 
— 500 
100 


|74-in -<liam = air-dielectric 
[ stub. Zo=50.0 ohms 


adjustable 


130 \3\s-in.-diam air-dielectric perforated 


70 J line. 7Zo=50.9 ohms 
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Ficure 10. Measurement results 
f equation (15) for connecting tr 
B, transformer 7 in tandem with trar 
with transformer 1 


unsformers listed in table 1: A, trans 


sformer 1; and C, transformer 8 in 


nd. The magnitude of correction that can be 
roduced by a given transformer may be estimated 
by considering, for example, transformer 8; a7=0.981, 
h’ 250 micromhos. Using eq (17) and (18), it 
ean be shown that a measured impedance of 
(00 —710.0 ohms at the input end of this connecting 
transformer corresponds to an actual load impedance 
of 291—j31.5 ohms at its output end. Thus the 
error in the resistive component is 3 percent, and the 
reactive component is in error by 3 to 1. 

The accuracy with which the values of a? and ab’ 
can be determined depends first on the absence of 
curvature in the plot of eq (15) as represented in 
figure 10. With no observable curvature, the pre- 
cision with which the straight line can be drawn 
depends on the scatter of the plotted points. The 
scatter, in turn depends on the accuracy with which 
the s and ¢ values are measured. The accuracy of 
the transformation corrections was +0.5 percent or 
better for a’, and +(5 percent+0.01 micromicro- 
farad) for ab’. 


5. Measurement of the Quality of Slotted 
Lines 


The n-s method may advantageously be applied 
to determine the uniformity of a slotted line and 
thus indicate the order of its utility, either as a 
measuring line or as a variable reactor. The aux- 
iliary equipment required in performing this quality 
measurement is an adjustable stub or a short-cir- 
cuited line stretcher connected to the output end of 
the slotted line. This stub need not be of high ac- 
curacy; its purpose is only to move a voltage node 
along the slotted line. The procedure follows the 
general n-s technique. At a half wavelength from 
input and output ends of slotted line, voltage nodes 
are selected to represent, respectively, So and 7». 
[t is desirable that the measurement frequency be 
high enough so that most of the line lies between S, 
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and 7). The movable short circuit is then used to 
place a voltage node at selected S locations, and 
the corresponding values of s and ¢ are accurately 
measured. The data are then plotted and the values 
of a? and ab’ calculated, by the usual procedure. 

These resulting corrections represent the total 
transformation produced by all discontinuities be- 
tween Sy, and 7). Usually it is found that the mag- 
nitude of these corrections is less than would be de- 
duced from the nonuniformity of the probe output 
voltage, indicated as the probe is moved along the 
perfectly terminated line. The choice of whether 
to use the n-s technique or the voltage uniformity 
method to determine the quality of a slotted line 
depends on the information desired. Variations in 
probe coupling do not necessarily affect the accuracy 
with which nodal positions can be located, but do 
result in probe output irregularities along the line. 
Hence, the n-s method (depending only on accurate 
nodal measurements) determines rather accurately 
the electrical uniformity of a slotted line, and voltage 
uniformity measurements (depending on line quality 
and probe coupling) determine the combined effect 
of probe penetration variations and electrical non- 
uniformity and establish the error limitations of 
complete impedance measurements. Thus, the n-s 
method offers means of separating the two major 
causes of slotted-line errors. It follows that 
slotted lines, which are poor as impedance-measuring 
instruments, often may serve satisfactorily as vari- 
able reactors. 

A way in which the uniformity of a slotted line 
can be qualitatively evaluated is by using it as a 
variable reactor in the n-s measurement of some 
discontinuity placed outside this line. Fixed 
continuities in the line, which change the electrical 
lengths of the s values, will cause curvature in the 
plotted results; when no curvature is obtained, the 
line is free from serious discontinuities. 


also 
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6. Conclusions 


The work described above was limited to lossless 
connectors. Transmission-line sections may also be 
considered connectors as long as_ their are 
negligible. A long lossless line may, of course, be 
looked upon as consisting of several connecting trans- 
formers in tandem; one can then obtain the trans- 
formation corrections of representative sections and 
calculate the transformations likely to result from 
the long line. 

The measurements were made on connectors 
having ideal junctions with uniform rigid-line sec- 
tions; therefore, some of the transformation cor- 
rections cited are probably better than those for 
connectors joined to commercial cables having 
braided outer conductors. One may, of course, 
measure a number of the latter samples and get 
comparative data for various types and makes of 
connectors, cables, and adapters, as well as for any 
combinations of these. The transformation cor- 
rections are obtained for discrete frequencies; it is 
pointed out, however, that in cases of rather small 


losses 





distributed discontinuities these corrections do not 
vary sharply with frequency. ‘Thus, measurements 
at several well-selected frequencies may be sufficient 
to cover a wide frequency range. 
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ontinuous Measurement of Atmospheric Ozone by 
Automatic Photoelectric Method’ 


Ralph Stair, Thomas C. Bagg, and Russell G. Johnston 


An automatic photoelectric instrument and method for the continuous measurement 
of the ozone in the earth’s atmosphere at low altitudes are described. The method is physical 
rather than chemical in character and is based upon the optical absorption characteristics 
of ozone in the Hartley and Huggins ultraviolet bands. The instrument makes use of a low- 
pressure mercury arc, which is situated at a distance of 1,450 feet from the recording station 
that employs a 1P28 photomultiplier as a detector. The light beam is modulated, at 510 
cycles per second, so that the output of the photomultiplier is fed into a tuned alternating- 
current amplifier and amplified to the recorder level. By means of a Geneva mechanism 
which changes the glass filters, the radiant energy from the lamp is separated into bands 
primarily at wavelengths 253.7, 365.5, and 405.0 millimicrons. From the ratios of the de- 
flections for the different spectral regions it is possible to determine ozone concentration in 


the range from a few tenths of 1 part to many parts per 100 million 


1. Introduction 


Although ozone is an important constituent of the 

mosphere because of its protective effect in absorb- 

» short-wave solar radiation, it is highly deleterious 
© rubber and related organic products. Because 
wone is a highly oxidizing or catalytic agent, the 
presence of organic particles in the atmosphere, to- 

ther with the effects of long-wave solar radiation 
and increased temperature at lower altitudes, results 
n the destruction of most of the ozone molecules 
below about 15 to 18 km under normal conditions of 
atmospheric circulation [1].2 As a result of the loss 
of ozone at low altitudes and a decrease in total 
pressure with altitude, the partial pressure of ozone 
near the earth’s surface is only about one-tenth that 
of its maximum at an altitude of about 22 km. Rel- 
ative to the total number of molecules of air per unit 
volume, these figures become about 2 and 600 parts 
per 100 million, respectively. 

Because ozone of significant amount is normally 
produced only at high altitudes in the stratosphere 
by short-wave ultraviolet solar radiation through a 
process of photochemical dissociation of the oxygen 
molecule, followed by recombination, any of the gas 
found at the earth’s surface must have been trans- 
ported there by some means [12]. Diffusion is slow 
and hence for the lower altitudes may be neglected. 
General large-scale air-mass movements associated 
with eyelonic and anticyclonic conditions, together 
with local air circulation resulting from thunder- 
storms and the like, must, therefore, be responsible 
for any ozone transportation at a rate sufficient to 
ipset the static balance at the earth’s surface. In 
fact, variations in the ozone concentration associated 
with air-mass movements are well known [13]. 

The primary purpose of this investigation has been 
to find a practical physical method for use in the 
routine measurement of the concentration of ozone 
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in the atmosphere at low altitudes, in particular at 
the surface of the earth. It may be that through 
such studies valuable information on air-mass move- 
ment or other weather phenomena mav be obtained. 
An instrument for this purpose should prove espe- 
cially valuable in the correlation of the deterioration 
of rubber products with ozone changes. 

Rubber, when exposed to ozone, develops a type 
of cracking of the surface not produced by other 
elements or materials normally in the atmosphere [2] 
These cracks develop much more rapidly when the 
rubber is cured under strain [3] or when it is under 
stress resulting from a bending or elongation of the 
sample. It has been found that in outdoor exposure 
of rubber the amount of cracking varies seasonally 
and with locality. The resulting cracking is suffi- 
ciently reproducible that it may be used as a rough 
quantitative measure of the amount of ozone by 
exposing 2 sample of specially prepared rubber [4] 
to the atmosphere under the tension produced simply 
by bending the sample. Although the results are 
greatly affected by temperature [5] the method has 
been found useful in rough quantitative determina- 
tions of ozone. 

By employing the rubber-cracking method of 
ozone measurement, a survey at a number of cities 
in 1949 showed variations from a severity factor of 
360 in Los Angeles and Denver to less than 30 in 
New York City and Conshohocken, Pa. [6]. These 
tests were in general agreement with rubber-deteri- 
oration experience for the localities involved. Nor 
is the local ozone concentration entirely dependent 
upon transportation of the gas from the strato- 
sphere. Tests made in connection with the Los 
Angeles smog [7] indicate the formation of excessive 
amounts of ozone by photochemical action associated 
with the oxidization of hydrocarbons by nitrogen 
oxides present within the area. 

In recent years interest in the deteriorating effect 
of ozone on rubber materials has increased greatly. 
The development of new synthetic elastomers having 
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varving degrees of resistant properties to ozone also 
has contributed to an increased interest in this field. 
Although ozone is only one of several agents influ- 
encing the deterioration of rubber materials under 
the complex conditions of exposure to weather, tt 
frequently proves to be the most severe [8]. Simple 
methods for the accurate measurement of | 
centration are urgently needed 

Various chemical methods for the measurement of 
ozone have been suggested and set up. Most of 
them are elaborate, require attention, and 
integrate over a relatively long interval of time. A 
method recently developed by Regener [1] has 
shortened the interval to about 30 min, as well as 
making all the processes of the measurement auto- 
matic. In the past, physical of 
ozone have been considered, but no real progress has 
been made in this direction, although the high optical 
absorption of ozone between 2,500 and 2,600 A lends 
itself ideally to the use of a physical (optical) instru- 
ment for this purpose. The present report deals 
with the development of such an instrument 
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measurements 
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URE | Optical density of ozone, 
Tsi-Ze 
[10] 


from the data 


rences 


and Choong Shin-Piaw given in refe 


2. Instruments and Methods 


The method chosen for investigation, as indicated 
above, is physical in character and is based on the 
optical absorption characteristics of ozone in the 
spectrum. Although ozone has several 
various wavelengths 


about 
ozone 


it does not extend to wavelengths longer than 
3,400 A, even for a thickness of 1 em of 
(normal temperature and pressure) which value is 


ultraviolet 
about 3 to 5 times the total in the earth’s atmosphere 


intense absorption bands at 


within the ultraviolet, visible, and infrared regions 
of the spectrum, the strongest is the Hartley band 


centered between 2,500 and 2,600 A. The optical 
density of ozone (designated absorption coefficient, 
a, [9, 10]) as obtained by Ny Tsi-Ze and Choong 
Shin-Piaw is depicted as a function of wavelength 
in figure 1. Although the absorption band is broad, 


at any time. 

As ozone has a maximum of absorption in the 
spectral region of 2,500 to 2,600 A a source having 
high optical emission at 2,537 A, such as a low- 
pressure mercury arc lamp, lends itself ideally to use 
in a physical setup. No other available light source 
has such ideal characteristics for this purpose 


TaBLe 1. Method employed in the calculation of ozone values 


The spectral transmittances of this ammount 


This example includes a single set of data for an ozone value of 1 part per 100 million in the path distance of 1,450 feet. 
and 


of ozone given in column 4 were calculated from the data of figure | and correspond to an equivalent layer of ozone 0.000442 cm (normal temp<ratur: 
pressure) in thickness} 
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the lamp is equipped with a glass tube that 
bs the 1,847-A emission line,’ very little radiant 
vy is emitted at wavelengths below 3,650 A 
r than at 2,537 A. The spectral-energy distri- 
on of such a source is given in table 1, 
or convenience of modulation by a sector disk 
a source concentrated in space is preferred 
ever, as po such source is commercially avail- 
a bank of five Westinghouse Sterilamps (non- 
1 type) was assembled, with a 2-in. section of 
1 tube (see fig. 2 for a block diagram of the 
p) as the souree of radiation? The radiant 
vy from the 2 in. these lamps is 
ulated by means of a rotating sector disk driven 
-hp, 1,800-rpm synchronous motor. The disk 


section ol 


\7 openings, so that the light beam is modulated 
One purpose of the modulated beam is 


10 es 
liminate as far as possible the effect of other 
ces of light that might fall upon the photo- 
etor. Other frequencies, in particular 120 ¢/s, 
e considered but were discarded in favor of the 
ev [11] 
The light source is electrically modulated at 120 
by the a-e line voltage. At first thought, it 
peared that this would be idea! for the purpose 
However, 60- and 120-cy electrical pickup and 
stray light modulated at the same frequencies ruled 
out the practical use of 120 c/s. The electrical 
modulation of the light source at 120 c/s has been 
found not to interfere with the chopped modulation 
at 510 c/s. Higher-frequency modulation would be 
preferred, except for the practical difficulties in- 
volved in the mechanical chopping of a beam of 
sufficient aperture to supply adequate radiant 
energy for the measurements 
At the detector station the modulated light beam 
is picked up by a 1P28 photomultiplier connected 
with a tuned amplifier and recorder. The special 
photometer unit contains a Geneva mechanism for 
automatically changing optical filters over the 
photomultiplier. A set of three Corning glass filters 
has been arranged on a disk, so that alternate 
readings are made at 2,537 A (principally), in the 
blue (including wavelengths 3,655, 4,050 and 4,358 
A), and at 3,655 A, (see fig. 3). Each filter is placed 
in front of the photomultiplier for 15 sec, after 
which the filter disk is rapidly moved to a new 
position for the next filter in the series. The fourth 
filler position carries an opaque shutter. Hence, 
during each minute a complete series of measure- 
ments is completed and recorded. 
By means of a timing clock the source is auto- 
matically turned on about 30 min before a similar 


' Since ozone is created from free oxygen atoms, which are in turn created by 

otodissociation in the Schumann-Runge continuum, starting around 1,900 A, 

the source will produce a large amount of ozone if the emission line 1,847 A is not 

wed. Such ozone production would greatly interfere with absorption 

easurements. In fact a small ozone concentration near the lamp results from 

’ A radiant energy [14]. This has negligible effects upon the measurements. 

No condensing lens or mirror was employed in this work. The light source 

allowed to shine normally upon the detector. Although the resulting 

sity was low, its constancy more than compensated for the reduced in- 

ty because of the refractive and reflective difficulties encountered in the use 

figured reflective mirror. As the source intensity and detector sensitivity 

not limiting difficulties in this work, precaution was taken to eliminate as 

y optical problems as possible. As a matter of fact, the size of the light 

e could be diminished appreciably without encountering any grave 
difficulties 
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“IGURE 3 Transmittance curves for the 
the relative spectral response of the 
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Corning filte re and 
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clock starts the receiving apparatus. I ikewise, the 
same clocks turn off the equipment at a prearranged 
time each night. As all the components are con- 
trolled by timing clocks, the equipment requires 
no personal attention during the periods when it is 
in operation. Only daytime servicing and noting 
pertinent weather data, etc., are required. 

The electric circuits for the 1P28 photomultiplier 
and for the tuned amplifier are displayed in figures 4 
and 5. As the radiant energy transmitted through 
1,450 ft from 2-in. sections of five Sterilamps is 
extremely low, the modulated source and the tuned- 
circuit amplifier were required in order to reduce the 
noise and to eliminate the zero drift resulting from 
the variable dark current in the photomultiplier. 
The photomultiplier was commonly operated at 
about 70 v per anode, which resulted in a voltage 
amplification of the detected signal of about 100,000. 
Following this, a tuned a-c amplifier (with a voltage 
gain of several thousand) was employed to improve 
the signal to noise ratio and to obtain sufficient 
voltage for operating a standard strip recorder. 

Because the photomultiplier is highly sensitive to 
changes of voltage, a carefully regulated power 
supply is essential for precision work. The power 
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Figure Electronte circuit employed in the photomultiplier unit, 
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Figure 5. Electronic circuit employed in the 510-cycle-per-second tuned amplifier. 
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ply used (see fig. 4) contains many of the circuit 
ments commonly employed in conventional units, 
keep the output voltage constant over a wide 
we of input voltage. As the electric load of the 
plifier is relatively constant, little attention was 
en to regulatory devices to cover load variations 
Che tuned amplifier (see fig. 5), peaked at 510 ¢/s 
certain characteristics that may be noted. The 
ut stage employs a 6J7 tube for producing a 
itively flat gain at a high signal-to-noise ratio. 
t little shaping is introduced by the coupling and 
dback elements in this circuit. The second stage, 
ploying a 6SH7 tube, contains a tuned secondary 
sisting of a variable inductor shunted by a fixed 
v-loss condenser. This combination results in a 
ss band of about The amplifier is stable 
all gain settings, vet has a sufficiently narrow pass 
nd to reduce noise and power-line frequency 
erference to a relatively low level. 
\mplifier gain is controlled by stepped resistor 
nks in each of the three stages. Ample shielding 
the first stage, together with plate decoupling 
iworks, reduces noise pickup and positive feedback 
a minimum. 
The output of the amplifier, after rectification by 
6H6 twin diode, is fed Into a commercial poten- 
ometer-type recorder For purposes of test and 
djustment, the a-c output of the amplifier may be 
read on a suitable vacuum-tube voltmeter. Thea-e 
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oltage output is linear for all values below about 20 
Recorded d-c values are slightly nonlinear as the 
esult of contact potential within the 6H6 diode 


ctifier However, this error becomes a second- 
order effect as ratios only are employed in the reduc- 
tion of the amount of ozone. If desired, any scale 
may be eliminated through an a-c voltage 
libration of the d-c recorder deflection. Usually 
lose linearity has been obtained for all deflections 
ibove 1 v when full recorder-scale deflection was set 

10-v (alternating current) measured at the 
est-meter position. 

The method employed in the evaluation of the 
amount of ozone in the 1,450-ft path between the 
ght source and the detecting system is illustrated 
by a sample calculation in table 1, which gives all 
the steps in the determination of the filter-trans- 
mission factors corresponding to a concentration of 

part of ozone in 100 million parts of air. Similar 
calculations were made for other concentrations of 
ozone. The results of these calculations applying 
to the particular photomultiplier and filters employed 
in the present work are illustrated by the curves of 
hgeure 6. 

The work is simplified greatly by the fact that 
relative values only of the lamp energy and photo- 
multiplier spectral response are required. Hence, 
amplifier and recorder sensitivities may be neglected, 
except to the extent that a reasonable deflection on 
the recorder chart should be obtained. 

It is recognized that low-pressure mercury-are 
lamps are affected in total and in relative spectral 
output by their operating temperatures. This is 
governed by the ambient temperature and especially 
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Fieure 6 Calibration curves for use in 


records into ozone concentrations in parts per 


molecules of atmos phe re 


by air currents flowing past the lamp tubes. Large 
changes in temperature of the lamp tubes were 
eliminated by compactly enclosing the five lamps in 
a box having a small quartz-glass window through 
which the radiation emerged. For best results, 
heater units, or even temperature-control mechan- 
isms, should be set up in the lamp box under situa- 
tions of extreme temperature variation, Otherwise, 
the temperature may be recorded and a correction 
made to the calibration factor 

Selective scattering of radiation from the light by 
the atmosphere has been neglected in this work 
Rayleigh scattering is small, and any error resulting 
from it is much than that caused by other 
erratic changes in the atmosphere. It 
that any scattering by dust particles is nonselective, 
and therefore does not affect the transmission ratios 
on which the evaluation of ozone is based. In order 
to correct for any scattered radiation from surround- 
ing city lights, the moon, or other sources, which 
might reach the photodetector, a shutter is auto- 
matically interposed in the light beam for about 
| min every 15 min. 

Absorption by various impurities in the atmos- 
phere, such as smoke or other combustion or chemical 


less 


is assumed 


vapors associated with industry or housing, produces 
errors of unknown magnitude. Hence, the method 
cannot safely be employed where any vapors having 
absorption in the ultraviolet spectrum (at the wave- 
lengths of the mercury emission lines) may be inter- 
mittenly present. Many chemical and combustion 
vapors are known to have high ultraviolet absorp- 
tion. Chimney vapors from the burning of coal, for 
example, when blown into the light beam cause 
erratic recorder tracings. 

The method is not suited to daylight operation 
without the installation of elaborate shielding as 
scattered daylight, principally in the spectral region 
of 3,200 to 4,000 A, will overload the 1P28 photo- 
multiplier and make its modulated response non- 
linear. When better filters are available, it may be 
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that the response can be confined to only those wave- 
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lengths that are not in sunlight. In that case, day- . 2.0 
time operation should be equally satisfactory. o aak 
\ | 
; 
3. Results of Ozone Measurements at i 
Washington, D. C. z 1.4 nr || 
z i.2b / \aprit 3, 1953 | \r 
Over a period of several months during the late = © \ | \ 
winter and early spring of 1953 an automatic record = 3 '°°f } | 
was kept of the ozone concentration over the grounds =z -8F / | 4 
of the National Bureau of Standards for a period of — 3 ¢} [ — ‘ 
about 3 hr during the early part of each evening. ° \ 
Figure 7 gives an example of sections of one of these “[ r } 7 
records made during an evening on which the amount & .2- N FEBRUARY 12,1953 \ 
of ozone varied from time to time. On most evenings ~ an are © Boge Re ee — Cate 
there was little or no change. Except for changes 8:00 8:30 9:00 9:30 10:00 o%3 
associated with incoming cold air-mass fronts, no EST, PM, WASHINGTON, D.C. 
ozone concentration exceeding a few tenths of 1 part 
in 100 million was normally observed. Associated 
with these fronts, a concentration of several parts Ficure 8. Variations in ozone concentration for two evenir 


1 , » . Ss ft) 1 oC a 
per 100 million often appe ared. Ometime s the con- February 12, 1953; mild, humid, cloudy, light northwest wind, follow 
ce ntration of ozone flue tuated (as show n in figs. ‘ passage of cold front at 1:30 p.m April 3, 1953; mild, sky clear, light south wi 
: preceding cold front and low moving in from northwest; possible upper 


and 8) or disappeared during the evening measure- _ instability 
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Fiaure 7. Sections of arecorder tracing for April 3, 1953, showing fluctuations in the ozone concentration at 
Washington, D. C. 
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ments. In all cases during the tests it had disap- 
eared by the following evening. 

Although the results are not conclusive, it appears 
hat the presence of a substantial amount of ozone 
t Washington, D. C. is of only temporary duration. 
fests with rubber samples substantiate this conclu- 
ion. No doubt the large amount of vegetation 
vithin the area is very effective in keeping ozone to 

neglible value except during the short time inter- 
vals in which cold air masses are bringing it to the 
earth’s surface within this area [15]. 

The high sensitivity of this equipment in the detec- 
tion of ozone suggests its use in precise routine 
measurements at various locations in any study of 
the association of ozone with rubber deterioration or 
with possible correlations between ozone distribution 

and other weather phenomena. 
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Applications of Dimensional Analysis to Spray-Nozzle 
Performance Data 


Montgomery R. Shafer and Harry L. Bovey 


Some possible applications of dimensional analysis in studies of the performance of 


continuous fuel-spray nozzles of the centrifugal tvpe are presented. 


Equations are developed 


showing the relations among nozzle capacity, mean drop diameter, spray angle, nozzle size, 


the density, viscosity, surface tension, and pressure of the fuel 


Using experimental data 


available at the National Bureau of Standards and in the literature, good correlation is 
shown in considerations of nozzle capacity, and fair correlation in those involving mean 


drop diameter and spray angle 


1. Introduction 


In the course of research sponsored by the Bureau 
of Aeronautics, Department of the Navy, on various 
fuel-handling and metering devices for aircraft 
engines, one of the accessories of interest has been 
the spray nozzle, which is widely used for the atomi- 
zation of liquid fuels in turbine-type engines. The 
performance of such engines is influenced by the 
characteristics of the nozzles, including the quantity 
of fuel delivered, the size of the droplets, and the 
angle of the spray. The latter, in turn, are known 
to be influenced by the configuration of the nozzle, 
the nozzle-pressure drop, and by certain physical 
properties of the fuel, including density, viscosity, 
and surface tension. 

The well-known advantages of dimensional analy- 
sis, namely, a reduction in the number of variables, 
the possibility that the effect of one quantity may be 
studied by varying another, and the ease with which 
it can be shown that certain variables are unimpor- 
tant, have not been utilized heretofore in published 
analyses of the performance of fuel nozzles. These 
advantages are demonstrated here by applying this 
method of analysis to the limited data available from 
research at the National Bureau of Standards and 
elsewhere on fixed, continuous-spray nozzles of the 
centrifugal type. 

Relations among several dimensionless variables 
are developed. The particular variables suggested 
herein are not new, and the reader will recognize in 
them the Reynolds and Weber numbers and products 
thereof. They seem convenient for showing both 
the individual and combined effects of nozzle size, 
fuel pressure, and the properties of the fuel upon the 
capacity of the nozzle, upon the angle of spray, and 
upon the mean drop diameter. 

The nozzles under consideration are of fixed con- 
figuration and deliver a continuous fuel spray. 
Atomization in such nozzles is produced by the con- 
version of potential energy of the high-pressured 
liquid to kinetic energy in the high-velocity liquid 
discharged from the nozzle orifice. Both the simplex 
type, having a single internal flow passage from the 
pressurized fuel supply, and the duplex type, having 
two separate fuel inlets and internal flow passages, 


are considered. Although only centrifugal spray 
nozzles are treated specifically, it is well known that 
the analysis holds for nozzles of other types and for 
restrictions through which a liquid discharges, sub- 
ject only to the limitation that the device have no 
moving parts. 


2. Nomenclature 


Dimensions 
(m, 1, t systems) 
D=a length characterizing the size of  (/) 
the nozzle 
d= mean drop diameter of the spray... (/). 
M = mass rate of flow . ((mt-!), 
AP= pressure drop through the nozzle (ml-'t-?), 
AP,=pressure drop through small pas- (m/-!¢-?). 
sage of duplex nozzle. 


-;=length ratios designating shape of Dimensionless. 


nozzle. 
mass density (ml-*). 
surface tension __. (mt~*) 


Dimensionless 


spray angle. 
(ml-"4t-!), 


absolute v ISCOSILyV 


3. Theory 


road 


A highly simplified description of the process of 
atomization from a simplex fuel nozzle, such as that 
shown schematically in figure 1, is adequate for 
present purposes. Pressurized liquid supplied to the 
nozzle inlet flows at high velocity through small tan- 
gential slots into the swirl chamber where a high 
angular velocity is attained. It then passes through 
the discharge orifice and emerges into the surround- 
ing atmosphere in the form of a hollow, conical sheet, 
the apex of which is at the nozzle orifice. The liquid 
sheet becomes progressively thinner as it moves away 
from the orifice. Eventually the sheet becomes un- 
stable, rupturing occurs, and small droplets are 
formed. 

The effects of the liquid pressure drop, density, vis- 
cosity, and surface tension upon the rate of discharge, 
the mean drop size, and the angle of the spray are to 
be proaes« Neat in this analysis. It is desired to de- 
velop relations by which experimantal data may be 
correlated to form characteristic curves and charts 
showing the relative influence of the various inde- 
pendent quantities upon nozzle performance. 
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Fiaurre | Schematic diagram of a wimpler spray norele. 


3.1. Flow Rate of the Simplex Nozzle 


Again referring to figure 1, the flow rate of the liq- 
uid through the nozzle could be influenced by the 
pressure drop, i. e., the difference in pressure between 
the nozzle inlet and the atmosphere into which the 
spray discharges, by the density, viscosity, and per- 
haps surface tension of the liquid, and by those 
linear dimensions defining the configuration and the 
orientation of the flow passages and swirl chamber 
Any one of these nozzle dimensions, such as the 
diameter of the orifice, may be designated as a char- 
acteristic length, J), and all others may then be ex- 
pressed as dimensionless ratios (7;,/,, r,) of this 
characteristic length. 

Assuming that all quantities that have an appreci- 
able effect upon the flow rate have been enumerated 
their intesedlntion can be expressed by the general 
equation 


F(M, AP, p.u.e,Dyry re, . © or) 90, (1) 


in which F merely denotes “function of.”’ 
By applying Buckingham’s pi theorem [1, 2], eq 
(1) reduces to 


1,=@(IM,, MMs,7,,72, . . «Pd, (2) 


in which each Il symbol represents a dimensionless 
variable or product formed through different combi- 
nations of not more than 4 of the 6 dimensional 


' Figures in brackets indicate the literature reference at the end of thi@ paper 
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quantities in eq (1). Also, each I is completely in 
dependent of the others, in that each contains on 
quantity not in either of the other two. 

The three dimensionless TI products may b 
formed by selecting any 3 of the 6 dimensional quay 
tities of eq (1) which are independent, and by associ 
ating these in turn with each of the remaining thr 
dimensional quantities. Since primary interest j 
in the effects of viscosity and surface tension upo: 
flow rate, it is desirable to derive each of the three TI 
products in such a way that one and only one will 
contain the quantities M/, wu, and « This leaves A?’ 
p, and )) as the three quantities that may be involved 
in all three products. On this basis, the quantities 
can be arranged in dimensionless groups such as 


Il, V/ 1h) Al’p; Il, Dy AP p/w; and Il DAP /o 

Substitution of these values in eq (2) gives 
M IP) Alp (1D, Al’p up, DAP/e, r,, Fr, . Pad 
3) 


The function is as vet unknown, and in fact it 
could have an infinite number of forms depending 
upon the shapes of the nozzles Howey er, if the con- 
sideration be now restricted to geometrically similar 
nozzles of fixed configuration but of different absolute 
sizes, i. e., to nozzles having the same r;, ratios, eq (3) 
reduces to 

WDPyAPp= (Dy APp/u, DAP/o), (4) 
in Which the form of the function ® will differ for 
various nozzle configurations 

As each of the quantities appearing in eq (4) can be 
measured, the explicit form of the function ® can be 
determined by experiment for any particular nozzle 
The experimental data might be plotted in terms of 
M/DP, APp as a function of DyAlp/u for different 
values of DAP/o. If such a plot shows that, regard- 
less of the value of the latter parameter, a single curve 
is defined, then it has been demonstrated that sur- 
face tension has little or no effect upon the flow rate, 
since o appears in eq (4) only in the product DA? /e 
Thus, to investigate the effects of changes in surface 
tension, it is necessary only to change the variable 
DAP/c. This can be accomplished by varying 
either o, AP, or D, although it is not simple to change 
D by fabricating geometrically similar spray nozzles 
of different absolute sizes. 

Herein lies one of the important advantages of 
dimensional analysis, in that the effect of variations 
in one quantity, surface tension in this illustration, 
may be investigated by changing some other quan- 
tity, in this case the pressure drop or the size 
Another advantage of this method is seen by com- 
paring eq (1) and (4). The former contains six inde- 
pendent dimensional variables, which have been 
combined into three independent products in eq (4) 
This reduction in the number of variables from 6 to 
3 simplifies the presentation and application of 
results. 
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3.2. Flow Rates of Duplex Nozzles 


As is shown in figure 2, duplex nozzles are pro- 


vided with two inlets and two separate internal flow 
passages. The inlet pressures may or may not be 
equal. From the viewpoint of dimensional analysis, 
the duplex nozzle of fixed configuration differs from 
the simplex nozzle in a single important respect, 
namely, that the pressure drops A? and A/’, for the 
main and small flow passages must both be included 
in the consideration Thus eq (1) for the duplex 


nozzle becomes 
F(M, AP, )” M, ¢, D, th, re,  v;) 0 (5) 


By forming the dimensionless groups as before, 
and by limiting the consideration to geometrically 
imilar nozzles, there results the relation 


V/ IP, Alp (1), Al’p M, AP AP,, DAP a) (t}) 


lt is known that variations in the parameter 
DAP /o have no appreciable influence upon the fune- 
tion ® for either the simplex or the duplex nozzles 
investigated. Thus DAP/e, and consequently sur- 
face tension, can be omitted from both eq (4) and 
(6 Equation (6) then becomes 


V/ 1P, Al’p (1D), AP p/u, AP/AP,). (7) 


In applying this relation to the experimental 
data, it will be found convenient to form charts of 
WibPyAlPp versus DyAlPp/u for various constant 
values of the pressure ratio A?’/AP, 
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Fiaure 2. Schematic diagram of a duplex spray nozzle. 


3.3. Mean Drop Size of the Spray 


Consider next the mean diameter of the droplets 
formed during the rupture of the conical sheet emerg- 
ing from the nozzle orifice. In an actual spray 
having a condensed volume, V, and consisting of NV 
droplets of various sizes, the mean drop size, d, is 
taken [3] as the diameter of a sphere having a volume 


Previous investigators [4, 5, 6, 7] have shown that 
the density, surface tension, and viscosity of the 
liquid, as well as the relative velocity between the 
liquid and the air just prior to atomization, may all 
influence drop size. Although the velocity of a 
liquid emerging from a centrifugal-type nozzle is 
not readily measurable, it depends upon the nozzle 
configuration and upon the flow rate, the latter 
being influenced by the pressure drop and by the 
density and viscosity of the liquid. The size of the 
drops also depends upon the shape of the spray 
chamber and the location of the drop-sampling 
apparatus within the spray. Thus, with the under- 
standing that the procedure applies only to geo- 
metrically similar nozzles, spray chambers, and 
drop-determination methods, the general equation 


F(AP, Pp, My, Oo, D,d)—0 (S) 


includes all the quantities that are expected to exert 
an appreciable influence when the spray discharges 
into still air 

Through the same procedure used in reducing 
eq (1) to eq (3), eq (S) may be reduced to the form 


d DD &(VAP oC, Doe ui). (9) 
Equation (9) is convenient for the present purpose, 
as if separates the effeet of fluid Viscosity If 


experimental results show that the effect of viscosity 
is appreciable, it would be convenient to form a 
chart of d/D versus DAP/o for different constant 
values of Dpo/u®. For any given value of J), the 
latter variable is constant for a given fluid at constant 
temperature, and can be varied conveniently by 
proper selection of the test fluids. 


3.4. Angle of the Spray 


Similar considerations may be applied to the spray 
angle. For this correlation, use is again made of 
the data of Rupe [3], who states arbitrarily that this 
angle shall be symmetrical about the axis of the 
spray and shall include 80 percent of the total- 
weight flow. It seems reasonable that this dimen- 
sionless angle, @, could be influenced by the con- 
figuration and size of the nozzle, and by the pressure 
drop, density, viscosity, and surface tension of the 
liquid. It does not seem probable that any other 
quantity will have a significant effect when the 
nozzle is discharging into stagnant air, except for 
the currents created by the spray itself. Thus, as 
before, the general equation may be written 


F(AP, p, uw, o, D, 0)=0, (10) 
which reduces to 
6= (DAP /a, Dpo/u*). (11) 
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4. Experimental Results and Discussion 


The validity of the relations developed above is 
demonstrated if plots of the experimental values of 
the dimensionless parameters give curves from which 
the deviations of the individual points do not exceed 
the experimental error. Such plots provide means of 
correlating the experimental data, and show the 
influence of the various quantities upon the flow 
rate, mean drop size, and spray angle. In_ the 
graphs that follow, the experimental data on the 
flow capacity of the nozzles were obtained at the 
Bureau. The data on drop size and spray angle are 
original with Rupe [3]. 


4.1. Nozzle Flow-Capacity Correlation 
The nozzles selected for the flow-capacity tests 
were types currently in use in aircraft gas turbines 
These included three simplex nozzles of one make, 
design, and nominal size, and two duplex nozzles of 
different designs. All were of the fixed configuration, 
continuous spray, centrifugal type, and each was an 
unmodified production unit 
Four different liquids were used in the flow- 
capacity measurements. The temperature of the 
liquid at the nozzle, as measured in each run by a 
thermocouple, varied from 77° to 90° F, depending 
upon the operating pressure and flow rate. The 
pertinent physical properties of these test liquids at 
the average operating temperature of 81.5° F are 
given in table 1 


rarie | Prope rlies of the capacity test liquids at S814 KF 
Description Den Viscos Surface 
sity. p ty. tension, ¢ 
ajicom p Dynesiom 
Commercial grade of heptane 0. 601 0.412 22. 8 
Light petroleum solvent Ts) nes 27.6 
45°, of white mineral oil plus 55°, of lieht 
petroleum solvent S13 2. 38 29.2 
7S. of white mineral oil plus 25°, of light 
petroleum solvent SS7 7.02 41.3 


The variations of the properties of the liquids with 
temperature over the operating range of 75° to 90° F 
were also known, and the actual values correspond- 
ing to the operating temperature observed in each 
run were used in computing the dimensionless 
variables. 

With a given liquid in the test system, the flow 
rates (\M/) were observed at nozzle pressure drops of 
20, 40, 70, 100, 175, and 275 Ib/in®. These flow 
rates were determined for each nozzle at each of the 
aforementioned pressure drops before changing the 
liquid. This procedure requires that each nozzle be 
installed and removed whenever the liquid is changed. 
Although the simplex nozzles were torqued to ap- 
proximately the same value at each insertion, it is 
possible that this process caused small changes in 
configuration, which would influence the correlation 
among results with different fluids. 
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The capacity of the duplex nozzles was determined 
in the same way as that of the simplex nozzles, ex- 
cept that the pressure ratio AP/AP, was controlled 
and measured. Ratios of infinity, 1.0, 0.6, and 0.0 
were selected arbitrarily for these tests. The 0.6 
pressure ratio had to be omitted in the case of nozzle 
2, because its flow rate proved so sensitive in this 
region that accurate measurements could not be 
made. 

In dealing with fuel nozzles it is common practice 
to express capacity in pounds per hour, pressure 
drop in pounds per square inch, density of the liquid 
in grams per cubic centimeter, its absolute viscosity 
in centipoises, and its surface tension in dynes per 
centimeter. Although these mixed units are used in 
this report for speciving experimental conditions and 
characteristics of the liquid, all measured quantities 
were converted to egs units before computing the 
dimensionless variables presented herein. Thus the 
dimensionless variables are based on flow rates in 
grams per second, pressure drops in dynes per square 
centimeter, densities of the fluid in grams per cubic 
centimeter, its absolute viscosity in poises, and its 
surface tension in dynes per centimeter. The dimen- 
sionless equations are independent of the system of 
units employed, and the user may select the system 
with which he is most familiar, provided only that 
the same units are used in developing numerical 
values of the parameters and in their subsequent 
applications. 
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In the present state of the manufacturing art it is 
not practicable to make the small and intricate flow 
passages of production nozzles geometrically similar. 
It was therefore anticipated that correlation of the 
data for the three simplex nozzles by a single curve 
would not be possible. Hence a separate charac- 
teristic curve is presented for each nozzle, and since 
D) for a given nozzle is constant, its value need not 
be determined. For simplicity the value of one 
centimeter is arbitrarily assigned to D. 

With this assigned value of D, the parameters 
MIP, APp and Dy APp/u can be evaluated from the 
experimental data for each nozzle and for each 
observed value of pressure drop, density and _vis- 
cosity. The results for the three simplex nozzles 
are shown in figure 3, and those for duplex nozzles 
1 and 2 are shown in figures 4 and 5, respectively. 
The same symbol is used in all these figures for a 
given test fluid. It will be seen that the results for 
each nozzle define a single curve or chart, from 
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which the individual points deviate by less than 

1 percent on the average. Such scattering ts 
commensurate with the experimental error, thus 
confirming eq (4) and (7 

Referring to figures 4 and 5, the capacity curves 
for the duplex nozzles show that the treatment of 
the two inlet pressures, through the formation of 
the dimensionless pressure ratio AP/AP,, is valid. 
Thus, through the relation expressed in eq (7), it Is 
possible to obtain characteristic capacity curves at 
selected, constant values of pressure ratio for the 
duplex nozzles. 

In general, the capacity curves show slightly 
increasing flow rates with increasing viscosity, for 
constant values of ), AP, and p. They also show 
a noticeable difference between nozzles that have 
the same nominal size and shape. After a curve of 
this type has been established for a given nozzle, it 
may be used for estimating the capacity of that 
nozzle to deliver other fluids. 

As an illustration of the utility of such capacity 
curves, the possible effects of changes in surface 
tension upon flow rate will be considered. In per- 
forming the capacity tests, the product DAP /o was 
varied between 4.5 10' and 83 10'. Despite the 
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fact that the parameter DA/’/o appears in eq (4), 
all values of 17/7, APp and of DAP p/u for a given 
nozzle lie on a single curve. As the surface tension 
does not appear in the plotted parameters and as 
the experimental values of the latter lie on a smooth 
curve, regardless of the value of DAP/o, the results 
indicate that the flow rate is not affected signifi- 
cantly by the surface tension of the liquid that is 
flowing. This conclusion is limited to variations in 
o that produce changes in VAP/o within the range 
specified above 


4.2. Mean Drop-Size Correlation 


Although no measurements of drop size or spray 
angle have been made at the Bureau, the applica- 
bility of eq (9) and (11) can be examined by referring 
to data obtained Rupe [3] presents 
extensive data on a few nozzles with different 
liquids, and his results will now be inserted in the 
equations. The pertinent properties of the liquids 
used by Rupe are listed in table 2 

Rupe’s method of determining drop size involved 
collecting the spray droplets for a predetermined time 
interval in cells having non-wetting glass bottoms 
These contained an immersion liquid having a lower 
density than that of the spraved fluid. The droplets 
collected on the bottom were photographed under 
magnification and their images were counted with 
an automatic electronic scanner-counter in 14 size 
groups down to 5 microns. The results of interest 


elsew here 


Tanie 2 Properties of the drep-size and spray-angle test 
liquids 
( poner I y» uaa ¢ 
Fluid N Spec Dope 
W ater Niger Add St : ray \ as 
tiv ! n ma] y 
nr ¢ Ww 
D 
wt ot ot 
1 “ma a7 oo “7.0 Lon ie) hg 
a a a7 0.70 is 1.0 1, At2 2;.8 
; vm 37 wi 7l alt 1. (48 {Kies 17.2 
‘ Loo a7 pAmN! 04 1. O88 whl u 
“18 a7 1.™ “0 1 om 1288 a2 
' “. 3 4.37 100 0 1a 1.284 21 
wo. eM i» a0 0. OOF 1. 420 13.6 
5 10D-octane gasoline v0 712 1. 407 M6 
t 
® Additive to fluids 2, 3, and 4 was glycerol; to fluids 6, and 7 was butyl 
alcohol 


Value of D taken as 1 cm in the parameter Dpe 
* Specific gravity and viscosity are estimated values 


here, as presented in figures 12, 13, and 15 of refer- 
ence [3], apply to a single spray nozzle with droplets 
collected from the same relative position in successive 
runs 

Rupe’s figure 12 applies for fluid 1, and shows 
mean drop diameter as a function of pressure drop. 
In his figure 13 the ratio of the mean drop diameters 
of fluids 5, 6, and 7 to that of fluid 1 are plotted as 
functions of surface tension at pressure drops of 
25, 30, and 100 Ib/in®. In his figure 15 the ratio of 








the mean drop diameters of fluids 2, 3, and 4 to tha 


of fluid 1 are shown as functions of viscosity a! 
pressure drops of 50, 100, and 150 Ib/in’*, ” Hi 
method of analyzing the results required data on tw: 
groups of fluids, one (fluids 1, 5, 6, and 7) havin, 
equal viscosities but different surface tensions and 
one (fluids 1, 2, 3, and 4) having equal surfac: 
tensions but different viscosities. Three separat: 
graphs were required for presenting the results, and 
these concealed rather than revealed possible une 
valuated influences. Had he planned to treat the 
results by dimensional analysis, there would hav: 
been much greater leeway in the choice of test 
fluids and in the selection of experimental pressure 
drops. The results could also have been presented 
on a single chart. 

In the present treatment the data of reference 
[3] have been substituted in eq (9), using the arbi- 
trary value of 1 cm for D. The resulting values of 
d/D and DAP /o are plotted in figure 6 for the seven 
test fluids. It will be noted that the points define a 
single curve, from which the deviations are usually 
within +10 percent. For values of DA?/o between 
2.810' and 1010', even better correlation is 
obtained. Considering the difficulties encountered 
in collecting, photographing and counting the 
numerous droplets ranging in diameter from 5 to 200 
microns, such correlation is surprisingly 
Nevertheless, an examination of the probable causes 
of the scattering is of interest, as it indicates the 
unevaluated influence in the 


good 


presence of some 
experiments. 

It seems improbable that the effects of viscosity are 
a primary cause of the scattering. Equation (9 
indicates that a plot of d/D) versus DAP/e should 
yield a chart with each curve representing a constant 
value of the product Dpo/y*. Thus all fluids having 
the same value of this parameter should define the 
sume curve, even though their separate values of 
p, o, and yw may differ considerably. Furthermore, 
if the plot of d/D versus DAP /e vields a single curve 
for any particular range of the product Dpo/y’, then 
it is known that changes in p or pw, which cause the 
product to vary within this range will not affect the 
d/D ratio, and consequently the drop size. 

Referring to figure 6 and table 2, the points 
representing fluids 2, 3, and 4 all define the same 
values of 











curve, even though their respective 
nN om oe ¢ am 
(Dpo/p?)X10-* are 27.8, 17.2, and 9.9. Thus, con- 
120 
FLUID NO, 2 3 a 6 7 
Vv | 
© 
2 
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40 
° 40 
(0 aps) 10-* 
Ficure 6. Mean drop-size curve (Rupe’s nozzle 4). 
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sidering only these three fluids, it could be concluded 
that variations in viscosity that cause the product 
Dpo/u*)X10~* to vary within the range 10 to 28 
will not affect the drop size. Considering fluids 5, 
6. and 7, for which the respective values of 
(Dpo/u?)X10~* are 28.2, 23.1, and 13.6, it appears 
that these fluids define a different single curve. As 
these two groups of fluids embrace the same range of 
Doo yp, the separation between the curves defined by 
the two groups cannot be attributed to either 
viscosity or density effects. 

Some possible explanations of the scattering in 
figure 6 are that some influential quantity has been 
omitted from the theoretical consideration; that 
the nozzle configuration may have changed, either 
from dirt or rust, during the experiments; or that 
properties of the test fluids at the instant of drop 
formation may have been different from the tabu- 
lated static values. Rupe mentions the latter 
possibility in connection with the surface tension of 
the water-alcohol mixtures. The method of collect- 
ing and counting the smaller droplets may also 
have been inadequate, thus resulting in appreciable 
errors in the total number of drops and consequently 
in their mean diameter 

The curve of figure 6 may be used to predict with 
fair accuracy the mean drop diameter from geo- 
metrically similar nozzles of different sizes, with any 
liquid and pressure drop, provided only that the 
values of DAP/e and Dopo/y? are within the ranges 
investigated. As the ratio d/) is nearly constant at 
the higher values of DAP/o, the mean drop siz 
would be expected to increase in direct proportion 
to the size of the nozzle. This is not true, however, 
in the region of lower values of DAP /c. 


4.3. Spray-Angle Correlation 


Spray angle data for a single nozzle at different 
pressure drops within the range from 5 to 100 
lb/in® are given in figure 24 of reference [3]. These 
data, obtained with liquids 1, 5, 6, 7, and 8 of table 
2, were presented in terms of @ as a function of AP, 
with a separate curve for each liquid. 

From these data, the dimensionless variables of 
eq (11) were computed, with the result shown in 
figure 7. Each symbol represents a single test liq- 
uid, and hence also a single value of the product 
Dpo/y®. As the individual points do not deviate 
from the curve by more than +5°, the correlation 
is considered reasonable. However, for DAP/¢ from 
3.510" and 10 10*, fluid 1 appears to define a 
different curve from that of fluids 5 through 8. As 
in the case of the mean drop-size correlation, viscos- 
ity effects are not believed responsible, because table 
2 shows that fluid 1 has a value of DAPo/,? within 
the range encompassed by fluids 5 through 8. 
Therefore, it appears that either the nozzle config- 
uration changed or some unknown influence has 
been omitted from the consideration. 
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Ficure 7 Spray-angle curve (Rupe’s nozzle 4 


5. Conclusion 


Dimensional analysis provides a useful and con- 
venient method for correlating and presenting exper- 
imental data on the performance of fuel spray noz- 
zles. It is shown that experimental values of flow 
capacity with liquids of various physical properties 
define a single curve for each nozzle, when plotted 
in terms of logical dimensionless variables. Cor- 
relations of available data on mean drop size and on 
spray angle are less exact than those on capacity. 
This lack of correlation may be ascribed, at least in 
part, to experimental error. However, the possibil- 
ity that some important quantity has been omitted 
in formulating the dimensionless relations is not 
excluded. 

The curves presented to indicate nozzle perform- 
ance apply specifically only to the particular nozzles 
for which they were determined. They are used 
solely to illustrate possible applications of dimen- 
sional analysis to the study of spray nozzles for such 
purposes as determining effects of individual vari- 
ables that cannot be changed readily and independ- 
ently; developing test programs leading to a desired 
end result from a limited number of experiments; 
analyzing, correlating, and interpreting experimental 
data; and predicting differences in performance due 
to the use of liquids of different physical properties. 
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Silver-Uranium System ' 
R. W. Buzzard, D. P. Fickle, and J. J. Park 


rhe phase diagram of the system silver-uranium was constructed from data obtained 


by thermal analysis, metallographic examination, and X-ray diffraction 
characterized by a eutectic having a composition of approximately 
and a monotectic occurring at 
The solid solubility of uranium in silver appeared to be between 0.1 


uranium occurring at 950° C, 


percent of silver 


The system is 
5 weight percent of 
1,132° C. and 0.23 weight 


and 0.4 weight percent at the eutectic temperature; no appreciable solubility of silver in 


uranium was noted, The temperatures of the 
tion of uranium apparently were unaffected by 


1. Introduction 


The purpose of this investigation was to make a 
survey of the silver-uranium system and to develop 
the constitution diagram by correlation of thermal, 
X-ray, and MIcroscopie studies 


2. Previous Work 


A survey of the literature showed that information 
on the constitution of the silver-uranium system was 
meager [1]2 It has been reported that silver and 
uranium separated into two distinet lavers on melt- 
ing, with very little evidence of alloving [2] and were 
mutually insoluble. A eutectic was observed in the 
silver laver, and considerable loss of silver by evap- 
orization was noted [2 


3. Preparation of Alloys 


The base metals consisted of ~99.9 percent of 
uranium (Mallinckrodt Biscuit) and 999-fine silver 
The allovs were prepared in beryllia crucibles in a 
high-frequency induction furnace under an atmos- 
phere of purified argon [3] 

It was observed from sectioned ingots prepared in 
the composition range of 30 weight percent of silver 
that silver and uranium separated into lavers when 
molten (fig. 1). Chemical analysis of a number of 
ingots revealed that the silver laver consistently 
contained 4 to 5 weight percent of uranium. This 
factor was later used to advantage by taking the 
4 to 5 weight percent of uranium alloy as a master 
allov; other alloys were made by direct alloving. 
The melting characteristics of the system required 
that the temperature of charges consisting of primary 
uranium and silver be raised above the melting point 
of uranium in order to obtain alloving. The molten 
silver had a high vapor pressure and, at the melting 
temperature of uranium, a high loss of silver by 
vaporization Was encountered. 

Uranium was an effective deoxidizer of silver so 
that, as uranium oxide was rejected from the melt, 
this oxidation caused a loss of uranium from the 
silver-rich alloys on each subsequent melting or heat- 





Investigation sponsored by the Atomic Energy Commission at the National 
Bureau of Standards 
? Figures in brackets indicate the literature references at the end of this paper 


gamma-beta and the beta-alpha transforma- 
silver 


treating operation. Analyses of the silver-rich alloys 
were made only for uranium and, as uranium oxide 
was included in the total uranium reported by this 
method, there is a possibility that sufficient uranium 
oxide was retained in the low-uranium alloys (<0.5 
weight percent) to cause the results of the chemical 
analyses for uranium to be higher than the actual 
uranium content of the alloys. In view of the change 
of composition of the allovs caused by these factors 
during subsequent thermal studies, the reported com- 
positions of the allovs (tables | and 2) were deter- 
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mined after examination rather than on the ‘“as- 
cast”? material, 
Tarte 1. Summary of thermal data obtained for silver 
iranium alloys 
Thermal arrest 
} Silver* | 
- 
| Liquidus | 64+1e—8 9 841 7 r] pa 
P) 
wt; oy ‘ ( ‘ ‘ | 
100 sei | . 
wow ane ais . 
01. ON un } ) 
nu @S7 une | . 
wah any aM | 
ws we WSS d 
ow st) ww . 
as. 70 way 770 Hs . 
o.28 1, 18 MA2 7 wl - 
(4 1, 128 any TAS ot * 
Os 1, 148 we Ti ow | ¥ 
0 1, 148 7#H2 an 
* Compositions show the ratio of the component metals and do not reflect 
the impurities 





Tari t 2 NSummar y of microscopic data oblaine d for silver- 
uranium alloys 
} 
Number of phases observed quenching temperatures, °¢ 
Urani . 
um * 
As cast aon 700 mM) wow om OSS | 
i 
- | 
0. 1 1 l ! 1 l 1 | 
a3 1 1 1 ! 1 1 1 | 
l 1 l l 1 ! ! 1 | 
‘ 2 2 2 2 2 Fusion Fusion 
43, 2 2 2 2 4 lhe 
18 2 2 2 2 2 Fusion le } 
2.1 2 2 2 2 2 do De 
2.0 2 2 2 2 2 do Deo 
0 2 2 2 2 2 do Do 
5.0 2 2 2 2 2 do Deo | 


* Compositions show the ratio of the 
the impurities 
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Fiaure | Ingot containing 30 weight percent of silver showing 


the uranium (dark) and silver (light) layers 


4. Experimental Procedures 
4.1. Thermal Analysis 


The thermal arrests were determined by time- 
temperature analysis, using a molybdenum-wound 
resistance furnace. The necessory apparatus in- 
cluded a program controller to maintain uniform 
heating and cooling rates and an electronic recorder 
to reproduce automatically the heating and cooling 
curves [4]. The thermal curves were obtained from 
specimens of approximately 50 g, which were heated 
and cooled at a controlled rate of 2° C/min with the 
furnace operating under an atmosphere of purified 
argon. The thermal arrests were derived from the 
cooling curves. 


4.2. Microscopic Analysis 


The specimens for microscopic examination wer 
mounted in Bakelite, ground on a series of silicon 
carbide papers, with the finishing paper of 000 grit 
and finished either on a levigated alumina lap o1 
electrolytically, using 70 v at 1% amp, 3 to 6 immer 
sions of 10-see duration in an electrolyte of the follow 
ing composition: 20 ¢ of chromic acid; 30 ml of 
water, and 100 ml of glacial acetic acid 

In general, the structures of the silver-rich alloys 
were developed by swabbing with a modified chromie- 
sulfuric-acid solution [5] composed of 1 part of a 
concentrated stock solution (100-ml saturated solu- 
tion of potassium dichromate in water, 2-ml satu- 
rated solution of sodium chloride in water, and 10-ml 
of concentrated sulfuric acid) in 9 parts of water 
The silver-rich alloys developed a grayish coating 
and required a light polish on a alumina finishing 
lap after etching. The structures of the uranium 
alloys were developed electrolytically in a 10-percent 
chromic-acid solution, using 0.5 amp and 4 vy at 
25° C 

Specimens. to be quenched were sealed under 
argon in high-silica glass tubes of 7-mm bore, homo- 
genized for 8 days at 925° C, and furnace-cooled 
The alloys were quenched in ice water, using the 
nickel-block technique, in the manner described for 
the titanium-uranium alloys {3}. 


4.3. X-ray Analysis 


The specimens used for the microscopic studies 
were subsequently placed in an X-ray spectrometer 
and a chart obtained of the (Cu-Ka) X-ray diffrac- 
tion lines at room temperature. By this method it 
was possible to identify the phases present in the 
alloys and correlate the X-ray and microscopic data. 


5. Results 
5.1. Thermal Data 


The results of thermal analysis (table 1) were 
were obtained from the two terminal-alloy series 
(0 to 0.23 weight percent of silver and 0 to 5 weight 
percent of uranium). The alloys in the composition 
range 0.23 to 95 weight percent of silver separated 
into two layers during the thermal-analysis cycle, 
each layer having an analysis within the respective 
terminal-alloy composition range, as indicated above 

Thermal-analysis data of the silver-rich alloys 
(1 to 5 weight percent of uranium) showed consistent 
arrests in the vicinity of 950° C. Similar arrests 
were observed for the uranium-rich alloys; in these 
alloys a second arrest was observed at approximately 
1,132° C. The uranium transformations were noted 
in all alloys and were apparently unaffected by 
silver. 

Two alloys (30 weight percent and 50 weight 
percent of silver) in the composition range of the 
liquid miscibility gap were prepared, and time- 
temperature cooling curves were determined on the 
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Arrests identical to those previously 
recorded for the uranium-rich terminal alloys were 
observed. The resultant ingots showed a = sharp 
line of demarcation between a silver-rich and an 


whole ingot . 


wanium-rich layer. The layers were parted and 
eexamined by thermal analysis. Each layer de- 
veloped arrests characteristic of the respective 
terminal-alloy series. The miscibility gap extended 
over the greater portion of the system. It was not 
possible to obtain liquidus determinations for com- 
positions in excess of 5 weight percent of uranium 
because the alloys separated into two layers during 
thermal analysis. ‘This indicated that the liquidus 
rose sharply from the eutectic to the monotectic 
temperature. On the basis of these observations, 
it was concluded that the composition range of the 
miscibility gap was determined to be approximately 
0.23 to 94.5 weight percent of silver at the mono- 
tectic temperatures. The vapor pressure of silver 
vreatly influenced the determination of the thermal 
arrests in the uranium-rich alloys. It was evident 
that silver had little effect on the thermal arrests 
of the uranium-rich alloys, and the temperatures 
used in the construction of the diagram were based 
on the values obtained for uranium. 


5.2. Microstructures of the Alloys 


In the study of the microstructures of these alloys 
table 2), using the number of identifiable phases 
present as the criterion, it was revealed that the 
solubility of uranium in silver lies between 0.1 and 
0.4 weight percent at room temperature, and little 
change tin this value was observed between 25° C and 
950° C 


as is “se Ud “Se 


Fiaure 2 





A eutectic-like structure was observed in both the 
uranium-rich (fig. 2A) and the silver-rich (fig. 2B) 
alloys. Results of the microscopic analyses showed 
the eutectic composition to be approximately 5 
weight percent of uranium and the monotectic 
composition to be approximately 0.23 weight per- 
cent of silver. In the composition range 0.4 to 5.0 
weight percent of uranium, some evidence of fusion 
was observed at the edges of the specimens quenched 
from 950° C and complete fusion was observed at 
955° C. This confirmed the eutectic temperature 
of 950° C as observed by thermal analysis. 


5.3. X-ray Identification 


X-ray diffraction charts were obtained, using the 
same surfaces of the specimens as were examined 
microscopically; the results of the X-ray and micro- 
scopic examination were then correlated. Only 
urantum and silver lines were observed for the alloys 
examined in this system, 


6. Proposed Diagram 


The thermal arrests, derived from the cooling 
curves, were obtained from “as-cast’’ alloys that 
were rapidly cooled from the liquid condition. 
Chemical analysis showed these ingots to be homo- 
geneous from top to bottom. Since both silver and 
uranium were lost during the thermal analysis cycle, 
the reported alloy compositions are those determined 
from the chemical analysis of samples cut from the 
ingot (at the thermocouple tip) after thermal 
analysis 
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As-cast alloys, showing characteristic eutecticlike structutes found in both the uranium-rich and the 


silver-rich terminal alloys 


A, 4.5 weight percent of uranium, chromic-sulfuric etch. B, 0.07 weight percent of silver, electrolytic etch, 10-percent chromic acid 


Magnification, X100. 
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Ficurr 3 Silver-uranium system, 94 lo 100 weight percent 


of silver 

Microscopic data one phase; ©, two phase; @, fusion 
Thermal arrests were observed at about 950 

762°, and 653° C in both the uranium and silver-rich 
alloys, and in the vicinity of 1,182° C in the uranium- 
rich alloys Microscopic studies revealed a eutectic- 
like structure in both series of alloys. It was indi 
cated that a monotectic occurred in the system at 
approximately 0.23 weight percent of silver and at 
the melting point of uranium; it Was not possible to 
differentiate between the melting point of uranium 
and the monotectic temperature. In the silver-rich 
terminal alloys, fusion was noted at about 950° C in 
the composition range 0.4 to 5.0 weight percent of 
uranium (table 2), indicating a eutectic occurring at 
950° C and approximately 5 weight percent of ura- 
nium, The microscopic studies of the quenched 
alloys used in the determination of the silver-ura- 
nium solvus were made on specimens that were heat 
treated prior to chemical analysis. It was well 
established that the solubility of uranium in silver 
is less than 0.4 weight percent at the eutectic temper- 
ature of 950° C (fig. 3). In evaluating the studies 
of the alloys in the composition range of less than 
0.4 weight percent of uranium, it must be borne in 
mind that the uranium-oxide content of the alloy 
appears as uranium in the chemical analysis, and it 
is quite possible that the solubility may be even lower 
than reported, No intermediate phases or inter 
metallic compounds were observed in the system 
as a result of X-ray, thermal, or microscopic analysis 
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Fiaure 4 Silver-uranium system 


7. Summary 


The silver-uranium system (fig. 4) was constructed 
from data obtained by thermal, microscopic and 
X-ray analyses 

The important features of the system include: (a 
A monotectic at 0.23 weight percent of silver and 
1,132° C with a liquid miscibility gap extending 
from 0.23 weight percent to approximately 94.5 
weight percent of silver. (b) A eutectic at approxi- 
mately 5 weight percent of uranium and 950° C 
(ec) The solubility of uranium in silver appeared to 
be between 0.1 and 0.4 weight percent at the eutecti« 
temperature while no appreciable solubility of silver 
in uranium was noted. (d) The gamma-beta and 
the beta alpha transformations of uranium were 
apparently unaffected by silver 


The authors express their appreciation to Martha 
S. Richmond for the chemical analysis. 
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Nonnegative Trigonometric Polynomials and 
Certain Rational Characteristic Functions 


Eugene Lukacs ' and Otto Szdsz ” 


Let 0b <b b, be n integers and 0< d,<d d,, be m(m~n) real 

numbers (not necessarily integers) Denote by g(@) the Vandermonde determinant formed 
m 

from the 67, bj, . . ., 62 with the first row replaced by 7 A, cos 6,0, where \ Ht (1 — b?/d? 
k=l 

for j}=1, 2, mn The question whether g(@) is nonnegative for all values of @ is closely 

connected with the problem whether certain rational functions are characteristic funetions 

Four configurations of the b,, , 6,, and d,, _ d,, are studied, which lead to non 


negative trigonometric polynomials 


1. Introduction 


In this paper we study several nonnegative trigonometric polynomials that are connected 
with the question whether certain rational functions are Fourier transforms of distribution 
functions (characteristic functions). 

We consider in the following a rational function ¢(f) whose poles and zeros have the same 
imaginary part @ and assume further that all poles and all zeros have the order one. Let 

id, Wy; i? ,(g— 1,2, . . . nv) be the poles and iW,. iw, (k=—1,.2, . . . sm), m<n, be the 
zeros of o(t), where v,;-a-+ib,, w,.—a-+-id,, and a >0, 0< b, < by: a Oe, OS CO ws. CO 
Such a function ¢(¢) fulfills the necessary conditions that a rational function must satisfy in 
order to be the Fourier transform of a distribution function. These conditions were given 
in two earlier papers {1, 2]. 

In the following we denote by 


6 hb h 
b?, b2, h 
V (1) 
g2e-2) fe-o. , ent 
the generalized Vandermonde determinant formed from the numbers 67, . . . , 62 with the 
exponents k,1,2,...,(n—1). We write further A, for the minor of the element in the first 


n 
row and jth column of this determinant, so that V,=>)(—1)’~'b?*A, 
l 


In order to decide whether the function ¢(f) is a characteristic function, we decompose it 
into partial fractions and apply the Fourier inversion formula. We see then after some com- 


putation that ¢(f) is the Fourier transform of a certain function f(z), 1. e., (0) | ef (x)dr, 


where f(z) is given by 
0, if r< 0 


i(r) 


A 


Ce “Ve ( 1)?-'A A cos by), ifx>0 (2) 


it 
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with 
n m 
a Il (a?+ 6?) Ud? 
0 j=l kel 
m n 
V, 1 (a? +d?) I 6? 
kel l 


» 


(2a) 
m }? , 
* nu 7B) (j=1,2,...n) 
If we introduce the polynomial 
Hl eC d; ( d2) a rd >I? - « (.3) 
where 
. | m . I 
‘ 1)" ¢ Id? nd ¢, B ot 
d, ( an { 1 
we can write 
n= FTG), 921,2,...,2 (4) 


The function f(z) as determined by (2) is real, after some computation it is seen that 


| f(x)dxr=1. If f(x) is nonnegative, it is the frequency function of an absolutely continuous 
0 


probability distribution. We obtain therefore the following condition: 
C'ONDITION l. The rational function o(t) is a characte ristic function if, and only if, 


(1 A, COs b,8), (1 A, cos 5,0), ae dA, COs 5,0) 
5? b?, petiecn b? 
g(@) > (5) 
’ 
be», ee bye 


for all @. 
If we consider a rational function without zeros, that is, if we assume that @o(¢) is the 


reciprocal of a polynomial whose roots are —ia, —iv, —i?, (j=1,2,..., ), then we 
derive by a similar reasoning the following statement: 


Conpirion Il. The function do(t) is a characteristic function if, and only if, 


1—cos 6,6, 1—cos 5,8, ..«., 1—cos 5,0 
5? b3, we tt b? 
Jo(9) , : eT - ‘ >0 (6) 
 aalinad rn ee oe G26 =! 


for all 0. 

This condition is analogous to (5), in fact, (6) may be obtained from (5) by putting A,=1 
for j=1,2, ...,n. We agree, therefore, to put 4,=1 whenever m=0. 

In this paper we consider the determinant (5) only for integer values of the 6, and show 
that it represents for certain configurations of the 6, and d, a nonnegative trigonometric poly- 


nomial. In this manner we can construct certain rational characteristic functions. 
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2. Some Lemmas 


For the discussion of the nonnegativity of g(@) we need the following lemmas: 
Lemma l. Jf b,, bo, . .. , 6, are integers, the determinant (5) can be factored so that 


g(@) (1 cos 0)"""A(ecos 6), (7) 


where A(x) is a polynomial in I of degree b, i--m, 


Li MMA 2. I} b,, b., owes b,, are odd entegers, and ifn m, we have 
q’ (0) =sin 0(1—cos* 0)"-""' B(cos 6), (S 
where B(x) is a polynomial of degree b, 2n-+-2m l 
LemMa 3. Let 
*1 
Px) (1—t?)"—' @* dt for s=0, 1, (9) 
Jt 
and 
Q(2,¢) =P, (x) —*Po(z). (10) 


The polynomial ()(x,f) vs nonnegative in the interval 1<x<1 7f, and only if, 


oS 2Z,. (11) 
Here z, is the root of the equation 
*1 
R(z)=Q (—2z,2) | (1—t*)""'(f— z*) dt=0, (12) 
which falls into the interval (0,1). 
l 


Moreover, let Wn \ n- > ome 
Then lim w,=p exists and is the root of the equation 
n- 
- 1. 
F(p)= 14 (3—p’) 20° | (1—®)e-"' dt=0 (13) 
J0 
which is located in the interval (0, y } ). 

The proofs of the lemmas 1 and 2 are completely analogous to the proofs given for the 
corresponding lemmas in an earlier paper [3], in which nonnegative trigonometric polynomials 
of the form (6) were studied. Lemma 3 is identical with lemma 3 of [3]; we omit therefore 
the proofs of the three lemmas and refer instead to [3]. 


The function g(@) as defined by (5) and (2a) depends on the parameters 


by, . . ., ba, dy, ..., dy. It is sometimes convenient to indicate this fact by writing 
g(0)=g(O\b,, . . ., Bardi, . . ., dm). Similarly, we shall write occasionally Vo= Vo(b,, . . ., 5.) 
for the Vandermonde determinant formed from the numbers 6?, . . ., 62. 


It is possible to obtain some of the results of our previous paper as particular cases of 
results obtained in the present paper. For this purpose we need the following lemma. 


Lemma 4. Let g(0\b,, . + b,:d;, . « « Gu) Oe th function defined by (5) and (2a). Tf 
we put d,=b,, we obtain 
q(6\ b,, e © 8 b,: di, ee = I> b,, dost; ° . ~~ 
Pattie < sen Wal 
a ee ee eee 14) 
Vo(d,, Oo, O41 b,) 


~~ PPPP SFP Peer ee) BP ke 


~ee SS © Bee ener & 














that is. the substitution d,—6, transforms the function 


oo eee eee 
VolOr. . s oo Oe) 


into a similar function that does not contain any more the parameters d, and 6, 


It is possible to give an elementary proof of (14) by a straightforward computation. 
The validity of (14) becomes apparent, however, if we remember that the rational function $(f) 
is the Fourier transform of Ce~“g(r\h,, eee ae . i) 


3. Nonnegative Trigonometric Polynomials 








In the following we assume that the 6,, d,, , 6, are positive integers, while the 
d,, .. ., d, may be any positive real numbers. The function g(@), defined by (5) is then 
a trigonometric polynomial. According to condition (1) the function ¢(f) is a characteristic 


function if, and only if, g(@) is nonnegative. The purpose of this paper is to consider con- 
figurations that lead to nonnegative trigonometric polynomials. For this discussion it is 


er 


convenient to introduce occasionally a new variable 


r= cos @, (15) 

and to use the notation | 
P(x) q(are cos 2), (16) 

Since we assume from now on that the 6,’s are integers, the substitution (15) transforms 


cos 6,6 into T,,(r), where 7, (x) is the 6,th Tehebycheff polynomial of the first kind. There- 
fore, we see from (5) that under the conditions of lemma 1, 





P(r) S35 (—1)? Ad 1 —Ay Te, (2) = A — 2)" A(z) (17) 
’ : 
In a similar manner, we have under the conditions of lemma 2, ‘ 
P"(z)e=S) (—1) Aa, Th (2) (1—z*)*-*"'B(z). (18) 
| 
When we make this substitution we shall use the fact that the trigonometric polynomial g(@) 
is nonnegative for all @ if, and only if, the polynomial (2) >0 for |r| <1. From (17) we obtain i 
unmediately 
Conpirion (111). The trigonometric polynomial g(@) is nonnegative if, and only if, the 
polynomial A(x), introduced by lemma 1, is nonnegative for |x| <1. 
We use condition (111) and the lemmas obtained in section 2 to derive conditions for certain 
trigonometric polynomials. The results are contained in the subsequent statements. d 


SrarementT |. Jf the 6, are the first n consecutive integers, i. e., by=j for j=1, 2, . . ., n, 
and if m=1, the trigonometric polynomial g(@) is nonnegative if, and only if, 2d} >n. 
In this case we have A, 1—/"/dj, j= 1, 2, . . ., » and obtain from (17) 


P(r) >» 1)" [1 (1- B) T2)| (1—zr)""“a+- br). 


From this equation it is possible to determine the coefficients a and 6. After some elementary 
computations we obtain 


[n(n I+ dil oy 
d? 


a nr? a? 


P ) on I . ) ; 
al - . )! n(n 1)+d? 75 
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Clearly, P(2)>0 for 2)<1 if, and only if, n?—d?)<n(n—1)4+d. This inequality can be 


written as —n?+n—di<n?—di<n’?—n+d}. Since n is a positive integer, this condition reduces 
to 2h? nh. 

COROLLARY TO STATEMENT 1. If the 5, are the first n consecutive even integers, i. e.. if 
b,=2) for j=1, 2, , n, and if m=1, then the trigonometric polynomials g(@) is nonnegative 
if, and only if, d?>2n. The corollary follows almost immediately from statement 1 


STATEMEN' 2 li the b are the first n consecutive integers, > @..-oe h } for } l, a ~ it, 
and if m=2, the trigonometric polynomial q(@) is nonnegative, if, and only if, a certain polynomial 
((r) is nonnegative for rr <l Tlere 

Q(2) = L(n,dy,ds)22 + M(nd,,d,) + N(n,d,ds), (19) 
where 
Li nd, dy) n' n'(d? t de) T Rd, / 
Vind, ds) 2n' On’ tn- n- n(d? ’ ke 2d, (20) 


Vind, d,) =n'—6n8 + 7n?—2n + (n?—n) (2+ 2) + Pe? 


PRooF From (2a) we see that 


forj7—1, 2, . . . n, and obtain from (17) 
"7 { F P\ ri) yn er 
P(x) 2 | 1) ‘Asp 1- (1- n)(1-a) TS (l—ar)" (ar? + br+e) 


From this equation it is possible to determine the coefficients a, 6, and c. After some elemen- 


tary computation, one obtains 


P(x) = “A (l—xr)" -*Q(2) 
did} : 
A particular case can be discussed easily. If d)=(n+-1) and d,=(n-+-2), we obtain from (20) 
L=4(n+1)(2n+1), M 4(2n4+-1)(2n?+-n+2), N=4n'+4n'+-19n?+-5n+4. It is then easv 
to show that Q(2) is positive in l<r<+1 so that P(z) is nonnegative for the 
values k,=j(j—1, 2, _,n) and m=2, d;=(n+1), d:=(n+2). This means in probabilistic 
language that the rational function with 2n+-1 poles —i(a+ij))(g—1, . . ., n) and —ia and the 


four zeros —i(a+id,), —i(a+id,) is a characteristic function. 


STATEMENT 3. lf the b, are the first n consecutive odd integers, — = 2) l 
for j= 1,2, .. . ,n, and if m=1 while n>1', then the trigonometric polynomial g(@) is non- 


negative if either one of the following two conditions is satisfied 


(1) d >y2n a 
(ii) d<y2n—1, 
but 


‘The case n= m1 is covered by Statement 1, and can therefore be excluded her 
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where 2, , is the root of the equation 


which falls into the interval (0,1 


If neither (i) nor (ii) is satisfied, then the function g(@) assumes also negative values 


COROLLARY TO STATEMENT 3. Jf for large n and d<_2n—1 


e 7) h) \(2n 1) d?| 2 


; » 
\ (2n—1)°—d? , 


O.59389157 ..., 


then q(@) is nonnegative. 
In formulating this statement, we wrote for d,; simply d. From (2a) we see that 
1)’. 
d l : for gam 1 BZ, ww cp 
d : 
According to lemma 2, as expressed by (18), we have 


> ' (d® — (27 —1)*)4 5725-1 (z) (1—2*)" B(x), 


where B(x) is a polynomial of the second degree. From this equation it is possible to determine 


\ 2n— I _(2n 1)*] 4) = 
Bia "yf d ] [1 d |: ¢? (21) 


Bix), we obtain 


where 


BO) [1 |, 0 ; (22) 
(2 )? 2 ¥ 
B(x) 27] 1 a z; B’(z) 2y -_ -. . 
d* rh J 


We first assume d>2n—1, then B’(x) <0 for O0<x<1. Since B(1)>0, we have B(x) >0 in the 
interval (0,1), from B(—z)= B(x) it is seen that B(x)>0 for |z|<1. Next we suppose 


(2n—1)>d>y2n—1, then B(z) has a minimum at the point r=0, so that B(x) > B(0O)>0. 





Summarizing, we see that B(z)>0 for |z|<1, provided d>y2n—1. But this means that 
P’(z) <0 in |z| <1. From (17) it is seen that P(1)=0 if n>m, so that P(z)>0 in |r| <1, and 
therefore g(@) >0 for all @, so that condition (i) is proved. 
We next consider the case d<y2n—1. We write 
2n—1)? (2n—1)—d 
B(x) l z?— ~ ‘ » 
th U (2n—1)?*—# 
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If we denote by 


and by 


then we have 
B(x) =y' (x? —% and P’(xr) ¥’(1—a2")""7(z 
Therefore, 


< 
P(r) v’| (1 —t*)""*(t y*)dt 


Condition (i) and the corollary to statement 3 follow then immediately from lemma 3; a brief 
table of the 2, is contained [3]. 

We proceed to discuss a more complicated case by assuming that the rational function 
¢(t) has two zeros (i. e., m=2) and that it has the first n consecutive odd integers as poles, so 
that 6,—2j)—1 for j=—1,2, ... ,n 

We obtain then from (18) 





P’ (x) y(1—2*)**(A+Cz?+ Dr‘), (23) 
where 
¥=2 2n—1)A, ) 
D-16,, 
1 (23a) 
Cc 1} (2n—9)Xr 2H L)A, il 
l= (2n L3n-+-22)r (2n—1)(2n—7)A,-1 n—1)(2n—1)A, ; 


We write ° 


( 
a » Ob A , 24) 
ir, 1X, 
so that 
P’ (zx) 4yX,,(1 —27)"-* (42 +- ax? —b (25) 


For n>m=2 we have P(1 0, and therefore ‘ 
* “1 n " 
P(x 4yX.,, (1—?t*)""*(4t' +-at?®?—b)dt forn>3 (26) 


Jt 
The nonnegativity of the polynomial ?(z) depends on the nature of the roots of the polynomial 
B(x) =42'+-axr*—b. This introduces the need to distinguish a number of cases, for this purpose 
we write A(r) =4r?+-ar—b=4(r—1,)(r—re), where 
T,==(—a—va’?+ 165)/8 
To ( a- ya t 166) Ss (27) 
B(x) =h(2*) 


In the following we denote by 


R,(2) | (1—t)* 20-9 (@— 2*)dt for 1=1,2. 


We are now in a position to formulate statement 4. 





5 We assume here A.0, We need not consider the case A, 0, since this reduces by lemma 4 to the case covered by statement 3 


* If nem, lemma 2 ts not valid, and moreover, P(1)>0. The case n= m<2 would therefore require a separate discussion 
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Srarement 4. Jf the 6, are the first n consecutive odd integers, that is, if 6,—2j—1 for 


ez). 2 _n, and if n>m=—2, then the trigonometric polynomial g(a) is nonnegative if, and 


only if, one of the following nine, mutually exclusive, conditions is satisfied 


(i) h, >O, a>0, and 6<0 

(ii) d, > O, a0, and a’? + 166 <0 

(iii) d, 0, a>0, 0< b<a+4, and Ry Vr.) — 1h (Vr) > 0 

(iv) 3 0, 0 >a 14+a>b aw/16, and Roly) mR (yr \>0O: if b< 0. we hare the 
additional condition Ry( — V1) —2hy(— yr) >0 

(Vv N O>a+S8>a+4>b a’ 16 

(vt) A, 0, 1>a 8, 0O>a4+4>b a7/16, and also Ry(Vr2)— mR (Vr) > 0; if 60, 
we have the additional condition Ry(— V1) — r2R\(— yr) >0 

(vil) \,<.0,a>0, and b>a+4 

(vint) r, <0, a0, and b> max |0, a+-4] 

(1X) 4, <0, a+ 4<b<0, and Ry(vr,) — nRy(vr,) <0 


The proof of this statement is rather tedious. It is, however, similar to the proof of state- 
ment 3 of the present paper and of statement (F) of [3] and is therefore omitted. 


We conclude this section by using lemma 4 to show that the earlier results may be obtained 
as particular cases of the statements of the present paper. We consider as an example 
statement | and write n=s+1 and put d,<k, where & is an integer. From lemma 4 we see 
that g(@\1, &+1:) and go(@\1, k—1, k+1, +1) are simultaneously non- 
negative so that statement 1 reduces to statement C of [3] 
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Flame-Emission Spectrum of Water Vapor 
in the 1.9-Micron Region "” 
W. S. Benedict,’ Arnold M. Bass, and Earle K. Plyler 


The emission of hydrogen-oxygen and acet ylene-oxygen flames has been studied by using 
a 15,000-lines-per-inch grating spectrometer with a lead-sulfide photoconductive detector 
rhe water-vapor band, which in room-temperature absorption extends from about 1.80 to 


1.96 microns, is broadened in emission to 1.7 to 2.2 microns 
measured, and a rotational analysis has been given for many of them 


Nearly 1,000 lines have been 
Many lines of the 


(011-000) band agree with those found in laboratory and atmospheric absoprtion, and many 


additional lines from higher energy levels may be identified 
(J +1)-u+y may be followed out to J=23 in both P and R branches 
lines in the upper-state bands (021-010) and (031-020) have also been identified 


The principal series J-y-——> 
Corresponding 
\ smaller 


number of lines may be assigned to the (041-030), (051-040), (012-001), 022-011), (032-021), 


111-100), and (121-110) bands, 


1. Introduction 


When hydrogen, hydrocarbons, or other hydrogen- 
containing compounds burn in air or oxygen, con- 
siderable energy is radiated in the infrared region. 
With spectrometers of low dispersion, three principal 
regions of strong emission were noted [1] near 2.7, 
1.9, and 1.4 microns. These correspond to regions 
of strong atmospheric absorption, due to water 
vapor [2, 3] and the emission must also be attributed 
to the H,O molecule. Studies of the emission with 
resolving power sufficient to separate lines due to 
individual rotational transitions in these rotation- 
vibration bands have not previously been reported. 
Due to the great increase in resolving power in the 
near infrared region resulting from the use of photo- 
conductive detectors, it is now possible to obtain 
spectra showing very much more detail, and to 
pursue the analysis to a point where the major part 
of the emission may be understood. In addition 
to the fundamental interest in determining the 
higher energy levels of the H.O molecule from 
spectra and in interpreting a rather complex spec- 
trum, the present studies are also of practical im- 
portance in providing the data from which relative 
intensities within the H,O emission spectrum may be 
used to determine whether the H,O molecule appears 
in combustion processes in nonequilibrium distribu- 
tions [4] and to determine the equilibrium “tempera- 
ture” of that molecule in flames. 

In this paper we present the results of measure- 
ments in the region 1.7 to 2.2 uw (5,900 to 4,600 em~') 
When acetylene-oxygen and hyvdrogen-oxygen flames 
were observed with a 15,000-lines-per-inch grating 
spectrometer in this region, the intensity was 
sufficient so that slits subtending a spectral interval 
as narrow as 0.2 em™' could be used. This region 
(referred to for brevity as the 1.9-~ region) is the 
most suitable for beginning the analysis of the 
vibration-rotation flame spectrum, for a number of 
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reasons. The sensitivity of the lead-sulfide detector 
and the efficiency of the spectrometer are high; 
at the same time, the intensity of the emission is high 
without involving emissivities so great that  self- 
absorption effects cause distortion of the relative 
intensities. (There is some distortion of the intensi- 
ties due to absorption of the lower-energy lines by 
the water-vapor molecules at room temperature in 
the spectrometer path, but this is confined to scattered 
frequencies in the center of the band, from 5,100 to 
5,500 em=' [5].) The vibrational transitions involved 
in this region are few and do not show strong per- 
turbations. The principal band, the only strong 
one in room-temperature absorption, is the transition 
(PP gh yh Pots) = (011-000), “Upper-state” — bands, 
e. g., (021-010), (031-020), ete., also are excited 
at the high temperatures of the flame; these also 
involve states whose energies are well known from 
absorption studies and cause relatively little con- 
fusion. The “convergence” resulting from the 
difference in rotational constants between the upper 
and lower states of all these transitions is such that 
the bands are well spread out, with no regions 
where many lines crowd together. 

Other regions in the emission spectrum of H,O 
have also been studied, but will not be reported in 
detail here. The 2.7-u region, in which the principal 
band is the », fundamental, is extremely intense, 
with the result that both self-absorption and atmos- 
pheric absorption cause difficulties. Moreover, the 
15,000 lines-per-inch grating with lead-sulfide de- 
tector loses sensitivity rapidly at wavelengths be- 
yond 2.6 uw. Hence only the short-wavelength por- 
tion of this band (2 branch) appears to advantage 
on our tracings. The 1.4-u region is quite strong, so 
that resolution only slightly inferior to that achieved 
in the 1.9 uw region was attained, but there are a num- 
ber of vibrational transitions in the region, and con- 
vergence is’such that the 2? branches are extremely 
crowded (forming a characteristic head near 7,500 
em”~'), while the P branches fall in the region of 
strong OH emission [6]. At wavelengths shorter 
than 1.33 yw, the H,O emission becomes progressively 
weaker. The present observations extend to 1 4, but 














slits as wide as | em™~' were required to develop the 


bands centering around 1.1 yg, so that little detail in 
the rotational structure could be separated. In the 
photographic infrared and red [7] good resolution 
has been attained, especially in the region 7,100 to 
8,800 A; rotational analysis of this interval is in 
progress, but the structure is more complicated than 
in the 1.9-u region because the vibrational levels are 
more numerous and show perturbations. In addi- 
tion to the emission measurements, studies have 
recently be made [8] of the absorption of heated 
HO in the long-wavelength region greater than 4 yu. 
The observations in the spectral regions other than 
1.9 w were useful in confirming many of the rotation- 
vibration assignments made in the present paper 


2. Experimental Procedure 


The spectrometer and the general arrangement for 


observing flame spectra have been described else- 
where [6]. In the present work, hydrogen-oxygen 
and acetylene-oxygen flames were observed, both 
just above the tip of the inner cone, and at greater 
heights in the outer cone. Various fuel-oxygen ra- 
tios in the neighborhood of stoichiometric were fed 
to the welding-torch burner. These varied condi- 
tions caused little change in the relative intensities 
of the emission in the 1.9-u region, a fact which in- 
dicates that practically all of the emission is due to 
the H,O molecule, and that the energy distribution 
for that molecule is not very sensitive to the burning 
conditions, This is to be expected, inasmuch as 
H,O is the major component of the burnt gas. Its 
emission, in regions other than those of active com- 
bustion (and with the type of burner used most of 
the radiation entering the spectrometer must come 
from such regions), would be predominantly thermal, 
characteristic of the weighted average temperature 
of the molecules in the line-of-sight of the spectrom- 
eter. With the fuel-oxygen ratios used, the maxi- 
mum (adiabatic) temperatures attained did not vary 
greatly, and were in the neighborhood of 3,000° K for 
both hydrogen and acetylene as fuels. The hydro- 
ren-oxygen emission was somewhat the stronger, the 
none proportion of H,O in the burnt gas presumably 
overcoming the slightly higher temperatures in the 
acetvlene-oxygen flame. 

Four of the hydrogen-oxygen tracings, showing 
optimum resolution, were selected for the measure- 
ment of wavelength. Wavelength standards were 
obtained by superimposing on the emission spectra 
lines of either neon or krypton from gaseous dis- 
charges, focused through the flame. <A smooth dis- 
persion curve was fitted to the standard lines, and 
the H,O lines were measured with reference to this. 
Prominent lines were measured on all four tracings, 
and were usually concordant to +0.2em~'. Weaker 
lines were measured on one or two tracings, using the 
strong H,O lines in the neighborhood as secondary 
standards. In the center of the band, the strong 
absorption lines due to atmospheric water vapor 
could be used as checks on the emission measure- 
Moreover, many weak lines that are ob- 


ments. 








served in the absorption of large columns of wat, 
vapor [5, 9] appear strongly in emission in the 
tracings because of the much higher temperature | 
the flame. The measured positions of these lin, 
are in good agreement with the absorption dats 
Because of the internal and external consistency « 
the measurements, it is believed that the frequenci 
of strong unblended lines are correct to within +0 
em~', and of the weaker or blended lines to +02 
em~', Wavelengths in air were converted to way 
numbers and corrected to vacuum in the usua 
manner. 


3. Results 


A complete listing of the lines between 4,600 and 
5,900 em~' has been prepared. Selected portions of 
this listing, corresponding to the regions of the P, Q 
and # branches as illustrated in figures 1 to 3, are 
presented in table 1. The first column gives the 
intensity, as observed on a typical H,-O, tracing, in 
arbitrary units of recorder deflection above the back 
ground of unresolved emission in the neighborhood 
The intensities are uncorrected for instrumental sen- 
sitivity and hence can be given only rough quantita- 
tive significance; the sensitivity, however, does not 
vary greatly over the frequency interval studied 
Lines whose intensity and position are very strongly 
affected by atmospheric absorption are marked with 
an asterisk (*). It should be mentioned that the 
peak deflection of the most intense lines is of the 
order of one-tenth of that for a black body at the 
flame temperature. Making allowance for the fact 
that the slit width exceeds the line width, it is esti- 
mated that the maximum emissivity at the center of 
strong lines does not exceed 0.2 to 0.4. This would 
appear to be experimentally confirmed by the fact 
that the peak deflection was very nearly doubled 
when a mirror reflected the flame upon itself, thus 
doubling the effective path of hot gas in the line of 
sight of the spectrometer. The second column of 
the table gives the position, inem~'. The remaining 
columns of the table list rotational assignments for 
the observed absorption. 

The measured region contains over 1,000 lines 
These include all lines with a measured intensity of 
10 or greater (the strongest observed line, 4,999.99 
em™~', having an intensity of 90), together with some 
weaker lines down to J=3 at the wings of the band. 
The noise level of the instrument and flame insta- 
bility corresponded to J(2, so that all lines, which 
have been observed on several tracings, are con- 
sidered definite peaks in the emission. Many weaker 
lines of probable but not certain reality might have 
been added, and there is also a background of un- 
resolved emission, whose intensity, in the stronger 
portions of the band, may reach 5 to 10 units. Some 
of the strongest atmospheric-absorption lines, arising 
from the lowest energy levels of the ground state, 
appear as dips in this background. There is little 
question but that if still higher resolving power 
could be used, many more “lines” would appear in 
this region. However, as will be discussed presently, 
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Fieurr 1 Emiasion spectrum of H.0O in the region 


the analysis has accounted for most of the expected 
transition intensity, and it is believed that a major 
fraction of the measured lines are single lines whose 


position and intensity will not be appreciably 
altered under higher resolution. 
Tracings of the observed emission, in three 


regions which best illustrate the nature of the 
spectrum and its analysis, are presented in figures 
| to 3. These are tracings of the spectrum to which 
have been added wave-number scales and various 
identifying marks pointing out some of the strongest 
The spectra show the degree of resolution, 
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and discrimination between strong and weak lines 
that has been achieved. Figure 1 covers the region 
5,055 to 4,908 em™~', which includes the strongest 
lines in the ? branch of the (011-000) and upper 
state bands. It is a region where atmospheri 
absorption is negligible. Figure 2, covering the 
region 5,450 to 5,280 em™~', shows how lines of the 
@ branch, which would normally be strong, are 
reduced in intensity by the atmospheric absorption 
The strongest dips in the base line are also due to 
atmospheric absorption, as may be verified by com- 
parison with the spectra in Nelson's report [5]. In 
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Figure 2 Emission spectrum of H,O in the region 5,280 to 5,450 em-' 
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Fiaure 3. Emission spectrum of H,O in the region 4,910 to 5,050 em 


figure 3, covering the region 5,680 to 5,510 em=', the microwave region (0.74 em™') to the green 
lines of the R branches of the (011-000) band appear visible (18,500 em~'), is well understood and quite 
essentially free from atmospheric absorption. These completely analyzed. The water-vapor molecule 
three figures include the region of strongest emission; — is an asymmetric-top rotator with effective reciprocal 
at frequencies greater than 5,680 em™' and less than moments of inertia, which in the ground vibrational 
1,910 em! the intensity falls off gradually. In the state have the values A=27.877 em™', B=14.512 
intermediate region 5,280 to 5,055 em™' the emission ecm™', C=9.285 em™', leading to an asymmetry 
is weak, since it lies near the band origins, and is parameter & 0.4377. It has three fundamental 
further reduced by atmospheric absorption. frequencies of vibration, two preserving the mole- 

cular symmetry » 3,657.05 em™', and » 1,594.59 


4. Discussion em~', and the third in which the motion of the 
H atoms is antisymmetric to the symmetry axis, 
The absorption spectrum of H,O, both in the »=3,755.79 em™'. In the absorption spectrum, 


region of 1.9 w now under consideration, and through- transitions occur from the — vibrational level 
out the rotation-vibration spectrum extending from to these three fundamentals, and to a very large 
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number of the possible overtones and combinations 
of them; in all, 43 bands have been observed. Many 
of these, however, are quite weak; in particular, 
transitions to levels in which ¢ is even are very 
much less intense than those in which ¢; is odd. For 
example, in the 1.9-~ region, there are the two 
bands (011-000), whose origin is at 5,331.19 em™', 
and (110-000), whose origin is at 5,235.0 em™'!: the 
latter is only about one thirtieth as strong as the 
former, and hence plays a very minor role in the 
absorption and cannot be detected in emission.’ 

In each band with Ar,-odd, the rotational 
structure is similar (except for differences resulting 
from the different values of A’-A’’, B’-B’, and 
c’-C”’ involved); namely, there are four strong 
branches, the relative intensity of whose lines may 
be caleulated according to the methods of Cross, 
Hainer, and King [10]. The Ry and Po branches 
extend, respectively, to higher and lower frequencies 
from the origin, and each consists of one line for 
each rotational level J) (—JsSrsJ). The Qs» 
and Qo branches fall near the origin (in the 1.9-4 
band, where A’-A”’ is positive, they lie mostly at 
frequencies higher than the origin), and together 


contribute one line for each J). (Q, for J—r even; 
Qi for J—r odd). The intensity of an emission 
line is given by the formula J=+'7, g e7*’""’, 


(in absorption, the formula is the same, except that 
v* becomes v, and £’ becomes /’’). ZL is the line 
strength, a number such that for each J,, 2L=2J’ 4 
1; for r~—J, L~J for Ry, and Po, and L ~0 for the 
@ lines: for r~+-J, L~O for R and P, and ~2. 
for QY, with a gradual transition between these 
limits for intermediate r. Since /” increases rapidly 
with r for a given J, the intensity of the R and P 
lines falls off rapidly as J++ increases. Since the 
Boltzmann factor and L oppose each other for Q 
lines, there are several lines of nearly equal intensity 
for each J, those with J—r low. Another term 
in the intensity formula is the nuclear spin factor 
g, Which is 3 for r odd, 1 for r even; hence the former 
lines are predominant. The final factor in the in- 
tensity is the vibrational transition probability 
T,, which in first approximation should be constant 
for all lines of a band, but which in some bands of 
H,O (but not 011-000) shows anomalies favoring 
one branch over another. 

The absorption spectrum in the 1.9- region is 
fully in accord with the above formula. It has 
been possible to identify all the lines down to those 
of quite low intensity, comprising the strong series 
Ry, Pa, Qn, and Qy» for BF’ s2,000 em. In 
addition, there have been observed lines of the 
series Ry, Py, and Qs, and @,; (these have much 
lower L but can be observed in the terrestrial at- 
mosphere, since the first two series lie at the extreme 
wings of the band) and, as was mentioned above, 
lines of the 110-000 band, which follow different 
selection rules, but which are all weak because they 
have a different and much lower 7,. 

* An exception to this statement holds for a few rotational levels of 110, which 


accidentally lie close in energy to levels of 011 of equal J and like total symmetry; 
transitions to these are strengthened by perturbation induced by Coriolis forces 





It is to be expected, and it has been found, that 
the emission spectrum differs from the absorptio: 
spectrum only as a result of the much higher valu 
of temperature in the Boltzmann factor term of th: 
intensity formula. 

The levels of higher energy that are sufficient), 
populated in the flame to emit with considerab), 
intensity manifest themselves in two different ways 
First, rotational levels of higher J} will be populated 
the fact that J increases along with J further aug 
menting the intensity of lines from such levels 
Hence, the (011-000) band will show many ney 
lines, especially in the long-wavelength (P branch 
and short-wavelength (R branch) wings of the band 
Second, higher vibrational levels, particularly th 
lowest ones involving », will be populated; sin 
7,, for transitions with Avy,=1, increases propo: 
tionally with », the intensity of all rotationa 
transitions in such bands as (021-010), (031-020 
ete., will be appreciable relative to (O11-000) at 
flame temperatures. These new “upper-state’’ 
bands will form sequences of decreasing intensity 
at longer wavelengths and hence will be mon 
prominent in the P branches than in the R, where 
many of their lines will fall in the region of atmos 
pheric absorption. 

A rough quantitative calculation of the expected 
relative intensities at flame temperatures has been 
made as a guide in the analysis of the new emission 
lines. For this purpose 7’ was chosen at 2,877° K 
(k7=2,000 E’). L values were obtained for J <12 
from the tables of Cross, Hainer, and King [10] 
interpolated to k& 0.437. For J>12, L was 
estimated by extrapolation of those tables, making 
use of the facts that the relative total intensities 
of the several branches are smooth functions of ./, 
that the distribution of intensity within each branch 
at any J as a function of r/J is roughly the same, 
and that the sum rules given in [10] must be obeyed 
The resulting L values may be in error by as much 
as 10 percent at J=20, but should be adequate 
for showing the trends of intensity. The energies 
for the Boltzmann factor were in part those ob 
served (see table 4), and in part were estimated 
for unobserved levels by extrapolation. The result- 
ing relative intensities, for lines of the Po branch, 
are presented in table 2. The intensity factors, 
gle-*"*", are given in a rectangular array, with 
J’ vertical and (J+r) horizontal. There are a 
number of cases, for high J+-r, and for low J++ 
at high J, where lines of odd and even r+ fall in 
unresolvable pairs; these are combined in table 2, 
and may be recognized, inasmuch as the 3:1 in- 
tensity alternation due to the nuclear-spin factor is 
absent. Similar calculations of intensities in the 
Ro, Qw, and Qo branches (of which the two latter 
also coalesce at high J+ 71) have also been made 
The calculated intensities, as listed for identified 
lines in the last column of table 2, are these rota- 
tional intensity factors for the case of the 011-000 
band (band a). For other vibrational transitions, 
an additional vibrational intensity factor must be 
included. 
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In table 2 a solid line has been drawn to delimit 
the range of J, values yielding Px intensities up to 
10 percent of that of the maximum. It is not to be 
expected that lines weaker than this could be identi- 
fied under the experimental conditions. A_ large 
proportion of the lines within this limit have been 
found, the exeeptions being principally among the 
levels of intermediate +r for which the energy does 
not vary with J and +r in as regular a manner as 
when r~—J or r~+J. The two other lines in 
table 2 at lower J similarly delimit the 10 percent of 
maximum-intensity range for two lower tempera- 
tures, 1,438° and 287° K, and illustrate how the in- 
crease in temperature results in the appearance of 
many new rotational transitions. 

The relative rotational intensities within the var- 
ious upper state bands will be nearly the same as for 
011-000, small differences arising from the different 
energies relative to the J=0 level. However, 
these may be ignored for the rough calculations 
made here. Table 3 gives the band origins and the 
relative intensities of the stronger expected bands, 
Loge = Tye (®- Bow /*T, where T, is the vibrational 
transition probability and (4,—,,) is the energy 
difference between the J=0 level of the upper state 
in the transition and J=0 of (011). The vibrational 


Relative intensities of type A P-branch emission lines of H,O at 2,877°K 


j f 
12 13 14 1 6 17 Is Ww 21-22 23-24 ma) 
iB 
t 
s Cc 
‘ ] 
’ 5 | 
\ 
! s | 1 
’ s 1 
s 2 l 
, ‘ j 
‘ 9 ‘ | 2 
l l i i } . 
l i l l ‘ 1 
i l ; l ! 2 
1 1 ! I 2 


energy levels involved, with a few exceptions, have 
been observed in absorption, and the corre sponding 
values of m=K—Et are given to 0.1 em"! 
unobserved levels are calculable from the vibrational 
energy formula [3], and should be accurate to a few 
wave numbers. 7, and J are expressed relative to 
1.00 for 011-000; the former is proportional to 775, 
a result that is theoretically exact to first approxima- 
tion, assuming that the combination transition (which 
is foridden if the dipole moment may be expanded 
in purely linear terms in harmonic normal coordi- 
nates) is due to cubic and quadratic anharmonic 
terms in the molecular potential function and to 
quadratic terms in the dipole expansion [11]. 

Each band in table 3 is designated in brief by a 
letter symbol. The calculated intensities are the 
products of rotational intensity factors like those in 
table 2 and the vibrational factors of table 8. A few 
of the higher bands in table 3 involve Darling- 
Dennison [12] resonance pairs (¢ 4202s With vP2r342). 
For these the intensities are taken proportional to 
the fraction of 2», in the perturbed wave function, a 
rather questionable assumption because the intensity 
relation in cases of resonance are complicated by 
phase factors [13]. Table 3 contains 11 bands with 
T>0.12. Table 1 shows that numerous assignments 


171 





have been made to the seven of these for which both 
upper and lower energy states had been observed in 
absorption; that the strongest lines of the P and R 
series have been identified for two bands, 022-011(j) 
and 032-021(k), of which the lower state energies 
had previously been known; and that a few assign- 
ments have tentatively been made for the strongest 
remaining band, 051-040(h), of which both the 
upper and lower state are newly observed. No 
assignments have been made for band i 


Tasie 3 Vibrational band origins and relative intensities at 
287? . 
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The identification of a line is considered sufficiently 
convincing for inclusion in table | when both the 
frequency and intensity checks are satisfactory. 
For a satisfactory frequency check, the line position 
must be within 0.4 em! of the calculated position. 
This calculated position is derived by taking the 
difference between the two energy states involved in 
the transition. The quantum numbers of these 
states are given in the central column of table 1. 
The corresponding energies are listed in table 4 for 
the two upper states of greatest importance in the 
1.9-u region, namely, 011 and 021. The lower-state 
energies of 000, 010, 020, 100, and 001 have been 
given in recent publications [8, 14, 15]. Other 
upper- and lower-state energy levels, particularly for 
the low J values, were derived from atmospheric- 
absorption data and will be published elsewhere. 
New levels of the upper states 012, 022, and 032 
derived from the present data are given in table 5. 
Agreement of an observed and calculated frequency 
implies that the frequency difference between two 
or more emission lines is in accord with the corre- 
sponding difference in two or more (often as many as 
50) absorption lines. It does not mean that the 
energy levels can be calculated to an accuracy of 
0.4 em” from molecular constants because the 
centrifugal-stretching and perturbation effects are so 
large as to make these calculations prohibitively 
difficult. However, the energy levels do agree 
roughly with what is to be expected, as may be seen 
from a comparison in table 6 of the rotational energy 
levels for J=12 of the 000 and 011 states, caleulated 





For 000 ¢)} 


on the rigid rotator approximation. 
constants used were A= 27.877, B=14.512, C—9.98 
em™'; for O11, A=29.54, B=14.59, C=9.014. Th 
calculated energies are compared with the observ; 
rotational energies E’—F, 9; the differences A a, 
chiefly due to centrifugal stretchings. They a: 
consistently greater by about 40 percent for 011, by 
follow a very similar pattern for the two vibration 
states. The only rotational level of 011 that show 
an anomalous Ais 12_,. The large discrepancy he: 
can readily be attributed to perturbation by 12_, « 
110, whose unperturbed energy should lie withi 
10 em™! of the unperturbed energy of 12_, of 01) 
because the calculated rotational-energy differen: 
is -+-94.3 em and the vibrational difference —96 
em”'. Similar calculations from inertial constan( 
have been made for other J values, and likewis 
confirm the assignments. 
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Another way in which the reasonableness of the 
observed energy levels may be seen (in addition to 
the good agreement of the combination differences 
and calculations of the type first mentioned) is in 
the arrangement of the energy levels, or of the 
observed lines, in series that show regular trends 
with J and (J+-r). Two such arrays will be demon- 
strated; others can readily be assembled from the 
dataoftable 1. In table 7 the lines of the Aa branch, 
for J28, are tabulated. Regularities may be noted 
for all the values of (J+ 7), and the way in which 
the various (J+-r) series diverge and converge is 
characteristic. There is an obvious perturbation of 
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the (.J+-r)=8 series at J=12, as mentioned above; 
similar perturbations of (J+-1)=6 and (J+1)=4 
occur at lower J values. The highest J entries in 
table 7 are rather uncertain, inasmuch as they can- 
not be confirmed by combination differences; on the 
basis of the series regularities and intensities they 
appear plausible however. Note that at frequencies 
greater than 5,560 em™' practically all of the identi 
fied lines belong to this branch, and that the array 
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of table 7 includes most of the strong lines in that 


region 

Table 8 shows the series regularities of the strong- 
est P= branches, those with J+7r=0, for all the 
bands; it will be noted that there are also regularities 
in the vibration sequences. In addition to the line 
frequencies, the line intensities are tabulated; when 
overlapping occurs, the intensities are given in 
parentheses. The theoretical rotational and vibra- 
tional intensity factors are also given 

Table 8 also demonstrates the second criterion, by 
which it is believed that the correct assignment has 
been made, namely, the agreement between the 
observed and calculated intensities. Both for the 
lines in table 8 and throughout the assignments there 
is gratifving parallelism between those quantities. 
Better agreement could hardly be expected, since it 
must be realized that the observed intensities are 
often irregualrly disturbed by weak neighboring 
lines (not to mention the large disturbances by atmos- 
pheric absorption of many lines, including all those 
of band a, between 5,100 and 5,600 em™~'), and that 
the analysis is incomplete, so that unidentified tran- 
sitions of higher J may be making additional con- 
tributions. The necessary condition for a satisfac- 
tory analysis, that no lines of high and medium 
intensity are predicted where there is no emission, 
has been met. 

The agreement between the observed intensities 
of unblended lines and those calculated for 2,877° K 
shows that there is a possibility of using the H,O 
emission for temperature determination. This cannot 
be done with complete confidence, however, until an 
even more complete analysis, or better resolution, 
has been achieved. On the basis of the present 
results it would appear that in the flames studied 





there is equilibrium within both the rotational ay 
vibrational degrees of freedom, and that the te: 


perature is within perhaps 500° of 2,877° K. Whi 
it is gratifying that an assignment can be given 
nearly 90 percent of the lines observed, and that th 
accounts for all the predicted lines of intensit 
greater than 10 and about 75 percent of the tot: 
predicted intensity, there are no clear-cut criteri; 
as there are in the spectra of diatomic molecule 
for determining whether a weak unassigned line may 
not be contributing to the observed intensity. Henc: 
although it is tempting to use the observed relatiy 
intensities of two lines that lie close in frequency by 
come from widely different energy levels, such » 
5,600.65 em~', (2’=8,920) and 5,602.55 em~! (B’ 
7,193) as a “thermometric pair’, it may lead to 
error unless still higher resolution reveals the ling 
to be truly single. 

The analysis also shows that there is no simpk 
correlation between the frequency and the uppe: 
state energy. In general, it is the case that the 
emission comes from higher and higher levels as on 
goes toward the wings of the band, but there ar 
frequent exceptions inasmuch as the various bands 
and the various sub-branches of a given band, over 
lap in a complicated manner. Hence it is question 
able whether measurements with low-resolution 
spectrometers, comparing the relative intensity and 
emissivity in frequency intervals 10 to 20 em™ wide 
could yield temperature determinations of high 
accuracy. There is no doubt however that measure 
ments of this kind would be useful in a rough way, 
particularly if the present results are used to choose 
frequency ranges where the radiation comes princi 
pally from levels of nearly the same /’ 

It was stated earlier that the empirical evidence 
indicated that the maximum emissivity did not 
exceed 0.2 to 0.4 at the cores of the strongest lines 
This is a condition that is necessary if self-absorption 
effects are not to vitiate intensity measurements and 
temperature determinations [16, 17]. The emissivity 
of the strongest line, at 4,999.9 em~', may be caleu 
lated approximately from the following data. The 
relative line intensity (sum of the three components 
listed in table 1) is 52. The total integrated intensit) 
of the band, as determined from absorption measure 
ments, is equivalent to 20 em~? atm! (atmosphere 
here refers to standard temperature and pressure 
The fraction of this in the 5,000 ecm line is 20 
52/0.0,, where the Q values are the rotational and 
vibrational partition functions, respectively 5,550 and 
2.62 at 2,877° K. Hence the integrated absorption 
coefficient of the line is 0.072cm™~?atm~'. If the path 
length is 2 em and the H,O content of the outer-con 
gas 80 percent, there are 0.80 273/2877 =0.152 
em atm of H,O in the line of sight. Hence the tota! 
line strength, S, is 0.072 0.152=—1.1 107° em 
The line shape will probably be Lorentzian, with a 
half-width y, of the order of 0.05 cem™. (The 
half-width of pure H,O vapor at 100° C and 1 atm 
is about 0.4 cm™'; its variation with temperature and 
pressure is uncertain, but the above estimate is of 
the correct order of magnitude, and is somewhat 
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vibrational intensity factor. > rif, rotational intensity factor 


greater than the Doppler half-width, y,=—0.023 
em”'). Henee the maximum emissivity, given by 
1—e *7,, would be 0.07, which is of the observed 
order of magnitude and definitely well below unity 
This rough agreement also tends to confirm that the 
excitation of the H,O vibrational states is predomi- 
nantly thermal In the 2.7-~4 fundamental, where 
the integrated intensity of the band is an order of 
magnitude greater, self-absorption effects would 
begin to be encountered. 
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